

REMARKS 



Applicants submit the following comments as a summary o£ the 
personal interview on October 3, 2001, and as a follow up to the 
Amendment filed Keener 5, 2001, in response to the outstanding 
matters. Additionally, the title of the invention has been 



amended 



Extension of Time 

A Petition has been filed under the provisions of 37 CFR 
SI. 136 for an extension of time to respond to the Examiner's 
Action of February H, 2002. The appropriate fee set forth in 37 



CFR 



§ 1.17 is filed herewith. 



Discussion of Space Occupied Rate and other features 

As requested by Examiner Sheinberg, the following comments 
are directed to a working example to explain the "spaced occupied 
rate" and other questions discussed with the Examiner. The 
working example is in line with claim 1. 

For the sake of convenience the sample molecule is a three 
atom model and is composed of three atoms 2a(Cl), 2b(C), 2c (C) as 
Cl-C=C. For the sake of convenience, the sample molecule is 
reduced into the two dimensional plane. 

(A) in items (1), (2), (3) and (4) of the attached 
sketches, the setting of the molecule surrounding surface, the 



r 





dividing of the molecule surrounding space into a plurality of 
component spaces, assuming of the frontier surrounding surface, 
and providing of the probe points on the frontier surrounding 

surface are shown. 

(B) Deriving of the space occupied rate (FFfield) is shown 

in the item (5) of the attached papers. 

The space occupied rate (FFfield) is given as the number of 
the probe points existing on the frontier spherical surface of 

each of the atoms 2a, 2b, 2c. 

Though both of the atoms 2b and 2c are in fact carbon atoms, 
it is assumed for the sake of convenience that all are probe 
atoms. Additionally, for the sake of convenience, all of the 
probe points are shown in the figure of item (5) . FFfield is for 
the atom 2a = 7, FFfield is for the atom 2b = 3 and FFfield is 
for the atom 2c = 8 are obtained. 

The fact that FFfield of the atom 2a or the atom 2c is much 
larger than FFfield of the atom 2b shows that the proportion 
occupied by the atom 2a or the atom 2c in the whole reaction 
characteristic of the molecule is much larger than that of the 
proportion occupied by the atom 2b. 

The fact that FFfield of the atom 2a and FFfield of the atom 
2c are not so different shows that the proportions occupied by 
the atom 2a and the atom 2c in the whole reaction characteristic 
of the molecule are almost the same. 
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e^r*>rc\-r (FFelectro) is 
■ i nn of the electrostatic factor iu*e 
(C) Deriving or tne cj. 

o£ the attached sketches, 
shown in the item (6) of (FFelect ro) , the 

In orde r to derive the electrostatic factor 

h 1 on energy between a unit charge set at the 
electrostatic interaction energy 

(rharae( 2a) ( charge (2b) and 
prob e point (i) and the charges (charge (2a) 

, a 2b 2C is obtained, and the 
charge (2c)) of the atoms 2a, 2b, 

v, opined is summed up all over the 
electrostatic energy thus obtained 

i-hP frontier surface of the atom 2a, and 
probe points on the frontier 

FFelectro of the atom 2a is obtained. 

2c is also obtained m the same 
FFelectro of the atoms 2b, 2c 

WaY ' „■ f.rtor (FF electro) of the one 

In fact, the electrostatic factor (FF 

, „, a ia large which shows that a 
atom, for example of the atom 2a, 

• , reacted with the present molecule, 

molecule, which is to be reacted 

oach to the atom 2a of the present molecule 
difficult to approach to the 

n has a positive charge and the molecule is easy 
the molecule has a posit 

o a if the molecule has a ney* 
t-n the atom 2a, n uuc 
approach to cne 

farihnr (FFsteric) in 
■ ■ ^-p t-Vipa steric tactor i^ au 
(D) Deriving of the swi 

r-n of the attached papers. 

■ the steric factor (FFsteric), the probe 
in order to derive the steric 

• set at the probe point (i) • The probe atom 
atom bp' carhon is set 

hnn is set on the probe point so as that the v 

with the probe 

radius of the probe atom sp' carbon is 

n .neraies (=steric interaction energies) 
point (i) ■ van der Waals energies 
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between the probe atom sp' carbon set at the probe point (i» and 
the atoms 2a, 2b, 20 are obtained by using a mathematical 
technique for calculating a van der Waals energy in a (well 
Known, molecular force field HM3 . The van der «aals energy thus 
obtained is summed up all over the probe points on the frontier 
surface of the atom 2a, and FFsteric of the atom 2a is obtained. 
FFsteric of the atom 2b, 2c is also obtained in the same 



way. 



The fact that the steric feature (FF steric) of the one 



atom, for example, of the atom 2a, is large shows that a 
molecule, which is to be reacted with the present molecule, is 
difficult to approach to (or 'attack to' ) the atom 2a of the 
present molecule since the steric hinderance due to the atom 2a 



is large. 



n n^raphprize the electronic and 
(E) In order to numerically characterize tne 

steric factors of the molecule, three Kinds of features (space 
occupied rate (FF field, , electrostatic feature (FF electro, and 
steric feature (FF steric,, are obtained for the atom 2a. The 
three Kinds of features are also obtained for the atoms 2b, 2c. 

Attached hereto as exhibits are the requested paper of a 
voronoi division. It is reiterated that the Voronoi division is 
not essential for the present invention, because the Voronoi 
division is only one embodiment of a procedure for spatially 
dividing the molecule surrounding space 21 into the plurality of 
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I 
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component spaces. For the Examiner's reference, Fig. 31 on page 
380 of the document shows similar procedure for dividing a space. 

It is respectfully submitted that the comments set forth 
above and those comments in the Amendment filed November 5, 2001 
fully respond to the questions raised during the interview on 
October 3, 2001. 



Conclusion 

in view of the foregoing, it is believed all the issues 
raised by the Examiner have been considered and appropriately 
addressed. It is believed this application is now in condition 
for allowance and action to that end is respectfully solicited. 

If the Examiner believes that a further conference would be 
of value in expediting the prosecution of this application, the 
Examiner is invited to telephone the undersigned to arrange for 
such a conference . 
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Attached hereto is a marked-up version of the changes made 
to the title by the current amendment. The attached page is 
captioned "v-~-inn wit h ^THnas to ^how changes made. " 



Respectfully submitted 
JACOB SON HOLMAN PLLC 



By_ 



400 Seventh Street, N.W. 
Washington, D.C. 20004-2201 
(202) 638-6666 
Date: April 15, 2002 
Atty. Docket: 5970/P63431US0 
JCH:DKD:VJB 




C . Hq>lman 
No. 22,769 
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Version with markings to show changes made. 

In the Title ; 

Please amend the title as follows: 

--METHOD FOR PREDICTING [FEATURES OF CHEMICAL] REACTION 
CHARACTERISTICS OF MOLECULES 
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m ! a suniev of a Fundamental Geometric 

if onoi Diagrams - A Survey 
SH? Structure 

|z AURENHAMMER ^ ScAteB ^ A.,, 

%e fur InformationsverarbeUung Techmsche 



This paper presentsasurveyo^^^^^^ 

rruXL in computational B"™^* St inside and outside computer ( science 
of the Voronoi diagram in a wide variety ot i particular emphasis on the 

°and surveys the history of SoSmic properties. Finally, the paper 

unified exposition of its ™ thel *£ c f *™ h g Qn Vo ronoi diagrams and related 
provides the first comprehensive bibhograpny 

StrUCtUreS ' ■ t a- F 2 2 [Analysis of Algorithms and Problem 

Categories and Subject ^Z^rr,- geometrical problems and 

Complexity]: Nonnumerica^ ^^^^ Combinatorics- combinatorial 
^ ion, ; W 1 V-Zr^Tc^on* Geometry and Object 

General Terms: Algorithms, Theory Marine combinatorial complexity, 

AdditionalKeyWordsandPhr^ 

convex hull, crystal structure, ^^^J^ irangement, fe-set, motion 
model , higher—^ proximity, random.ed 



Introduction 

Jgmputational geometry is ^err^ 
Pith the design and analysis of afco 
v^Ithms for geometrical P^^^d 



— ^Hp^s 0 ?computer sSence-such as com- 
^■?teuter graphics, computer-aided design, 
fcloboticf pattern recognition, and opera^ 
HlE- reseLh-give rise topro^njthat 
lifeherently are geometrical This is one 
^tgfeason computational geomet ry has at 
^flracted enormous research ^terertuiihe 
: feast decade and ^ a well-estabh^ar^a 

^Ittoday. (For standard sources, we reto to 
Hthe survey article by Lee and Preparata 



^Jaaers familiar with the literature £ 
computational ^^J^^Sb, an 

SUed the ^^^S 
trend can also be obser ™i i 
rial geometry and £ * journals 
ber of articles m natural sci ent«J 

Sat address the Vo^ JW™^ 




j^ine survey — . - 

permission to copy without fee all or ° f * t he ACM copyright notice and thattta* for £ ompu ting 

•fe distributed for ^^gi^uiSl copying is by P^^^Sm. 

JSfcnd its ^ appear, and no^e ^ gi .^ & fee and/or specie p 

£®^»Machinery. To copy ^erwiae,« ^WSOl 50 
il© 1991 ACM 0360-0300/91/0900-0345 S0l.au 

■^.■l'j! 

ACM Costing Survey, Vol. 23, No. 3, Septet 1991 
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developed for the computer constructicp 
and representation of Voronoi diagraipsp. 

The intention of this survey is threj^ 
fold- First, motivated by the fact thaM^ 
Voronoi diagrams have been (re)invente|g| 
and studied fairly independently m t|J||^ 
applied natural sciences, in mathemat^ 
ics, and in computer science it present^ 
sketches of their historical development 
in these three areas. Second, it surveys^ 
i-i. 4- — a nr. Voronoi diagrams andMSIl 



diagrams, 
Basic Properties of the Voronoi Diagram 



rule- Each point is associated with 
the region of the plane closest to it, 

^ny do Voronoi diagrams receive so 
much attention? What is fecial about 
Zs easily defined and visualized con- 
sSuct' lt seems three main reasons are 
struct, it Voronoi diagrams 

iSTfnaturfin various situations. In- 
deed several natural processes can be 
2 to define particular c asses of 
Voronoi diagrams. Human intuition . * 
Ljp- ffuided by visual perception. It one 

c ^underlying structure, the whole 
sees an unaenyms ° v,icrKpr 
situation may be understood at a higher 

tpresting and surprising 
pTopertils; for instance, they «™ ^ 
to many well-known W^^*^ 
tures This has led several authors to 
beheve that the Voronoi diagram is one 
of L most fundamental constructs de- 
fined by a discrete set of points. Finally, 
Voronoi diagrams have proved to be a 
pow°eXl too? in solving seeming y 

tS attention of computer scientists 
in the last few years. Efficient and 
easonably simple techniques have been 



tarv though important, properties, ui "«s«f. 

that will «^ 
feeling for this structure. We also mtr 0? ^ 
notation used throughout this pa^ 
5? He also Preparata and Shamos^ 
f5851or Edelsbrunner [1987] for source^ 
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of the 

set of n points icaueu ~ --- jgp 
For two distinct *. 

p as to q. Formally, 
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dom(p,<z)M-im*.P)^<)||| 




Computat 

Irith the 
lithms f oi 
lition, otr 
areas of c 
|uter gn 
"robotics, 
^ftions rese 
P&herentl 
^reason c( 



w S denoting the euclidean distanc|^^^ 

fiction Cleafly, dom(p,q) is a /°*^f 
SS Xne bounded by the perpendicular,^ 

nan pidnc v. v.icppfor sepavM^ 

from those closer to .9 anQ J™" rpffiori 0 f a«? ^ . - :«ast deca 
the separator of p and q. There^ K^W^ (F 

site peS is the portion ^^gsitehe surv. 
ing in all of the dominances c 
remaining sites in S. Forma 

reg(p)= 0, dom(p.q) 

qeS - { p) 
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Voronoi diagram for eight sites in the 



ce the regions are coming from m- 
-ting n - 1 half planes, they are 
x polygons. Thus the boundary of a 
a consists of at most n - 1 edges 
mal open straight-line segments 
ertices (their endpoints). Each point 
t edge is equidistant from exactly 
ites, and each vertex is equidistant 
at least three. As a consequence, 
.gions are edge to edge and vertex 
rtex, that is to say, they form a 
mal partition of the plane. This 
ion is called the Voronoi diagram 
of the finite point-set S (Figure 1) 
e that a region, say reg(p), cannot 
mty since it contains all points ot 
ane at least as close to p as to any 
3 ites in S. In particular, p e reg( p). 
ows that V(S) contains exactly n 
s Some of them are necessarily 
nded. They are defined by sites s ly- 
the boundary of the convex hull ot 
mse just for those sites there exist 
arbitrarily far away but still clos- 
»fo vertices occur if and only it all 
n S lie on a single straight line, 
legenerate configurations also im- 
! existence of regions with only one 
inded) edge, Otherwise, three or 
•dges meet at a common vertex. It 



iuex hull of S is the smallest convex poly- 
contains S. 



should be observed that each vertex is 
the center of a circle that passes through 
at least three sites but encloses no site 

Although n sites give rise to ( 2 ) - 
0 (n 2 ) separators, only linearly many 
separators contribute an edge to V(S) 
This can be seen by viewing a Voronoi 
diagram as a planar graph with n re- 
gions and minimum vertex degree 3^ 
lach of the e edges has two vertices, and 
each of the , vertices belong s toat 1 Bart 
three edges. Hence, 2c > 3w. Euler s re 
lation n + u - c > 2 now implies e S Z n 
- 6 and v < 2n - 4. Thus, for example, 
the average number of edges of a region 
Ses nora D chieve six; there are less than 
3° edges, and each of them belongs to 
exactly two of the n regions. 

The linear behavior of the size of the 
Voronoi diagram in the plane means that 
roughly speaking, this structure is not 
much more complex than the underlying 
configuration of sites. This is one of the 
ma n reasons for the frequent use of 
Voronoi diagrams. A second reason is 
that V(S) comprises the entire proximity 
formation about S in an explicit and 
computationally useful manner. For ex- 
ample, its applicability to the important 
^ office problem (see below) is based 
on the trivial observation that a point x 
falls into the region of a site p and 
only if P is closest to x among all sites in 
S Moreover, if site p is closest to site q, 
Then reg(p) and reg( ? ) share a common 
■ edee This particularly implies that tne 
closest pair of sites in S gives rise to 
some edge of V(S). 



Applications in Computer Science 

To substantiate the usefulness of the 
Voronoi diagram in computer science we 
briefly describe four situations where this 
structure is used. The practicality and 
diversity of these applications will 
impart the appeal of Voronoi diagrams. 

Associative File Searching 

Consider some file of n 

records referring, for example to latitude 

and longitude of a city or to age and 

ACM Computing Surveys, Vol. 23, No. 3, September 1991 
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income of a person. Suppose we .are ^given 
an additional target record R and we 
quickly want to retrieve that record of 
the file that matches R the best. If re- 
trieval occurs frequently on the same file, 
a supporting data structure is called tor. 
T^ssocLve file searching problem 
was first posed by Knuth [1973] a a 
two-dimensional generalization of the 
usual (one-attribute) file-searching prob- 
lem It possibly is best known in its geo- 
metric version under the name post-office 
problem: Given a set S of n sites m the 
plane (post offices), report a site closest 
to a given query point q (the location of a 
perso'n). Note that there exists a trmal 
anytime solution by computing all n 
distances. There are various f go"thms 
in computational geometry that need 
the post office problem as a subroutine 
[Preparata and Shamos 1985J. 

Shamos [1975a] first observed the 
relevance of Voronoi diagrams to this 
problem. A site p is closest to q if and 
onlv if q falls into the region of p. In a 
Reprocessing step the Voronoi diagram 
of S is computed. To report a site closest 
to a it now suffices to determine the 
region that contains q. For this so-called 
point-location problem, efficient solutions 
have been developed by Kirkpatnck 
fl983], Edahiro et al. [1984], and Edels- 
brunner et al. [1986]. In particular, point 
location in a Voronoi diagram with* 
regions is supported in 0(log n] > time ; and 
0(n) storage overhead. This shows that 
the post-office problem can be solved by 
means of Voronoi diagrams m logarith- 
mic query time and without increasing 
Se order of space, which is well known 
to be optimal already for usual file 
searching. 




sir* 



Figure 2. Voronoi diagram for two dense clusters. ^g^igf 



attribute data, similarity is reflected by 
the proximity of sites in the plane. Prox- 
imity, in turn, is revealed by properties 
of the Voronoi diagram for these sites. 
For instance, dense subsets of sites give 
rise to Voronoi regions of small area; 
regions of sites in a homogeneous cluster 
will have similar geometric shape; for 
clusters having orientation-sensitive 
density, the shapes of the regions will 
exhibit a corresponding direction sensi- 
tivity Figure 2 gives an example. Anuja 
[1982] showed how to use these proper- 
ties for clustering and matching sites. 

Voronoi diagrams support various 
clustering techniques used in practice. 
What is required at any stage of 
the clustering process is often little 
Ire than the retrieval of the nearest 
neighbor sites of specified sites. They 
can be reported easily by examining the 
edges of the regions of the specified sites 
This applies to several hierarchica 
methods [Murtagh 1983], partitions! 
methods [Asano et al. 1988], and meth; 
ods involving cluster selection [Aggarwai 
et al. 1989]. 
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Cluster Analysis 

The problem of automatically clustering 
data arises frequently [Hartigan 19751. 
Finding clusters means determining a 
partition of the given set of data into 
subsets whose in-class members are simi- 
lar and whose cross-class members are 
dissimilar according to a predefined simi- 
larity measure. In the case of two- 



Scheduling Record Accesses 

Consider a mass storage system repre- 
sented by a two-dimensional array of grid 
noints each capable of storing one record. 
The toe ^ the read/write head to move 
from point (x,y) to the point ■ 
proportional to \ x - u\ + \ y - v\, ™v 
distance between these points measured 
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Figu re 3. Voronoi diagram and minimum span- 
Zg trees in the L.-metric. 

in the ^-metric (i.e., Manhattan dis- 
~f One nroblem that arises m ac- 
batched requests is the following: 
TTi Ae head movement to be sched- 
S n order to retrieve the requested 

nq801 A satisfactory approximate solu- 
ofcan be obtained * ^£ 

voronoi diagrams. In an initial step the 

Voronoi diagram 

SP S e v An edge of the tree can only 
quickly. An eog regions are neigh- 

X Lee and Wong 1980] Figure , 3d- 
ustrates this diagram (solid) & and the 

SThtl is 
Sun a Set of 2 of "^sHu 
tree edge is traversed twice to access a 

the sites. 



Figure 4. Power diagram for seven circles. 



Collision Detection 

An important topic in controlling the rno- 
tion of robot systems is that ot c0 ^ sl0n 
detection. For a robot moving in ^ an ob- 
stacle environment, one needs to deter 
Se collisions between moving robot 
auhnarts and stationary obstacles or be 

wSlwo separately moving ; subpart of 
the robot. In particular able to 

turn is important in order to be able 



. , Wore a collision will 

stop the system betore <* 

have occurred. environment 
The robot system and the em 

on the *^^J&ZS££ 
ttK^^^^lue^ of the 
Sm T°h f^^g£ 

may be assigned to the circle 

^iLoTcircle^ ^ stemming from 

sections (among circu; 

the same [1985 , 

however if the po Unes 
Voronoi diagram ^ ^ is 

(Figure 4). 
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1. HISTORICAL PERSPECTIVE 

The history of Voronoi diagrams can be 
traced back to the middle of the nine- 
teenth century. Although the spectrum 
of scientific disciplines that include , in- 
terest in Voronoi diagrams is broad, three 
aspects have been emphasized: 

(1) Their use in modeling natural phe- 
nomena 

(2) The investigation of their mathemat- 
ical in particular, geometrical, com- 
binatorial, and stochastic properties 

(3) Their computer construction and 
representation 

Accordingly, Voronoi diagrams are 
useful in three respects: As a structure 
per se that makes explicit natural pro- 
cesses as an auxiliary structure for in- 
vestigating and calculating related 
mathematical objects, and as a data 
structure for algorithmic P™ bl ^ s .^ 
are inherently geometric. In all three 
applications, efficient and practical algo- 
rithms for computing Voronoi diagrams 
are required. Since the first application 
bears the initial seed for their investiga- 
tion let us consider the role that Voronoi 
diagVams play in the natural sciences 
first. 

1.1 Natural Scientist's Viewpoint 

To visualize the appearance of a Voronoi 
diagram in nature, one could think of the 
three-dimensional space being subdi- 
vided into a manifold of crystals: From 
several sites fixed in space, crystals start 
erowing at the same rate in all directions 
and 'without pushing apart but stopping 
growth as they come into contact. Lhe 
crystal emerging from each site in this 
process is the region of space closer to 
that site than to all others. In other 
words, the regions form a Voronoi 
diagram in three-space. 

1.1.1 Domains of Action 

Most of the early work on Voronoi dia- 
grams was motivated by crystallography. 




Figure 5. Cubic crystal structure. 



The objective in this respect was the 
study of regions arising from regularly 
placed sites. Figure 5 illustrates some of 
the regions defined by 15 cubically or- 
dered sites. In his comprehensive article 
on crystal structures, Niggh [1927] calls 
them Wirkungsbereiche {domains of ac- 
tion), a term widely used nowadays and 
used with many synonyms even m the 
1930s as demonstrated in the clarifying 
note by Nowacki [1933]. Much effort has 
been devoted to the fundamental crystal- 
loiaphical question: Which types of do- 
mSnS of action are capable of filling the 
plane or three space completely if only 
congruent copies (and certain motions 
are to be used? Significant work was 
done, among others %">*&^ i 
Delaunay (Delone) [1932], Nowacki 
fl 976^ by Koch [1975]. In fact, the 
question is a deep mathematical one, and 
Te will come back to it in the discussion 
of mathematical aspects of Voronoi 
diagrams. 



1.1.2 Wigner-Seitz Zones 
A physicochemical system consists of a 
number of distinct sites typically 
molecules, ions, or atoms^ Equilibrium 
and other properties of the system de- 
pend on the spatial distribution of the 
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„hich can be conveniently repre- 
dividing the space between 
Wording to the nearest-neighbor 
r/esulting Voronoi regions are 
called Wigner-Seitz zones after 
rand Seta [1933], who were the 
\ use them in metallurgy. Frank 
Lner [1958] use these zones m the 
ga P tTon of complex alloy structures, 
19701 bases his survey of cubic 
Structures on an intimate y ^ 
oncept. Allotting zones to sites is 
Jest to molecular physicists bio- 
tf material scientists, and physr- 
Sdste Via this approach Brostow 
SSe [1975] estimate th< , coardnuj 
amber of liquid argon, David and 
r 19821 study certain solvation 
j es Bamberger and Goodisman 
[nterpret the small-angle acatter- 
"Systs, and Augenbaum and 
[1985] treat large-scale hydro- 
ic codes-just to name a few. 

■ohnson-Mehi and Apollonius Model 



earest-neighbor rule (or eqmva- 
the crystal growth model men 
earlier) forces the Voronoi regions 
convex polyhedra. Allowing the 
s to start their growth at differ* 
gives rise to hyperbohcally shaped 
I The resulting Johnson-MeU 
was proposed by Johnson and Meh 
as a more realistic model of struc 
,r minerals. Figure 6 pvesan j- 
.ion; crystals are augmented with 
.ate of birth. , , 

rials growing ^^X/icaX 
srent rates give rise to sphencaUy 
I regions forming the Apollonius 
This structure can also be ob- 
as cell structures of plants or m 
ma de out of soap bubbles jMatzke 
Tst\ev 1946; Smith 1954; Williams 
It further appears as covering ar- 
' plants and as areas of best- 
3d transmitters [Sakamoto and 
i 19881- Weaire and Rmer [1984] 
comprehensive survey and various 
ical data. Interestingly, the : equi_ 
n state of a spider web constitutes 
>ralized, although still polygonal, 



j- o^, This follows from in- 
Voronoi diagram^ This toiio 

Since between ^"^teethey 
responds to the tension of the edge they 

define. 

7 1 4 Thiessen Polygons 

Graphics, interest £ Vg^-^ 

who proposed the name i 

Ins 'creating Thiessen wi* the.ide. 

VJTs ttyl* r esin geographical 
a p S attern P analysis, as o^«and?or 

T^ons o urLn planning), Gambini 
fS ] and Boofs [1W91 (market areas 
LiybftJ anu nolygons are used 

generalized Thiessen poiyg ; Mehl 

al. [1988] ^^^fj^al^g^L 
of Thiessen polygons from ajeogr P 
and economic point of view, see 
[1986], Eiselt and Pederzoh [1986], and 

Sib Suzu£ 97 and In [19861 report on the 
given subdivision of a geogr P ^ 
S^IXbSTS^ in the optimal 
JSW foW«- IIS 

l"^AnTare^ minimum spanning 
fr:f ^shortest connection network 
rS^re 8) the Delaunay triangulation 
fnd tL Gabriel graph. These three 
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Figure 6. Johnson-Mehl model. 
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concepts are intimately plated to 
Thiessen polygons as will be discussed 
later. 

1 1 5 Blum's (Medial Axis) Transform 
Yet another name for the Voronoi dia- 
gram is popular in the natural sciences. 
E's transform of a set of sites. ^Con- 
cerned with biological shape and visual 
science Blum [1967] used it for modeling 
new descriptors of shape. In general the 
main problem in pattern recognition is to 
Xact characterizing elements from a 
gSen (site) pattern. When no geometry 
5al model of the pattern is available its 
SruXral description may be based on 
il p motion of neighborhood of a site ana 
^ in particular, on Blum's transform^ 
F«SiSd [1979] successfully compared 



n Q731 curvature properties of a given 
£ correspond to topologica proper- 
ties of its transform (compare Figure »)• 
S smooth pieces of the contour (solid) 
SSiUbe viewed as sites of generali ed 
Sane By definition, each point of the 
tSorm (dashed) is equidistant from at |f 

least two sites. In this context, the trans- ^ m 
S i is called the rnedtal «w or m l,< 

skeleton of the contour. The medm ax* 

ESS /nd Lantuejoul and 
Maisonneuve [1984]. f annl : ra . 

Although the foregoing ^^al 
tions of Voronoi diagrams m the natural 
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r 



,. Minimum spanning tree for the nine 
litals of Austria. 

» is not complete, we refrain from 
'details and turn now to their role 
mathematical environment. 

thematician's Viewpoint 

ry earliest motivation for the study 
3 noi diagrams stems from the the- 
• quadratic forms. Gauss [1840] 
ed that quadratic forms have an 
retation in terms of Voronoi 
ms for sites that areparallelohe- 
ordered in space. His idea was 



_ » - - . «■ 



\ J 



Figure 9. 

1968]. 



Contour 



and its medial axis [Philbrick 



or S voromi diagram, and xndeed 

ACM Computing Surveys, Vol. 23. No. 3, September 1991 



if 

ar 



m 



■ m 

i as 



I 



, , .•.••.OHVi.SVJ 



ti>*t:v, r 



,;*v*<v So '.Y'v'- v.- » '■■ ■•'"V ■ 

1 S'-^i 1 -.'^' 1 ''.'' 'V.'. 1 V-vi* 








354 ' Franz Aurenhammer 
these terms are the most popular ones in 

^^^^ b r 

grams and their S enera . th undfir . 
SjjAfJSU Propels 
and applications), as well as points cm 
Sain features of Voronoi f ^ms that 
have potential applications but may be 
not familiar to computational geometry 
researchers. 



1.2.1 Regularly Placed Sites 

Mathematical interest in Voronoi dia- 
lviamema divided according to 

^Xr'egutly or irregularly placed 
whether regular y where 

sites are involved ^ main 

. a. filing of R is a cov- 

by cfosed :i\s whose inter, 

ors are pairwise dis 

by convex polyhedra with the propen,y 
tLt the group of motions mapping the 

interest. 1 heir poiy ster eohedra 

dLaS by Schonflies [189U One ^ of 

the central ^"J^^^^Ste 
tallography has been ™ hem 
all stereohedra and to clB ^*^ w 
according to their crystallography 

progress was made for the sub- 
clas of V plfswhedra. They can be mtar- 
preted as Voronoi regions and were called 
Leciol fundamental domains by 
ScSie/ [1891] and — J octoj 
bv Niggli [19271. Laves [19301 and 
D y elaunay et al. [1978] enumerated _ all 

Delaunay [1932] showed that the num 
ber of facets (faces of dimension d - D ot 





Figure 1 0. Plesiohedron with 18 facets [Grunbaum 
and Shephard 1980]. 



a stereohedron in K d is fmite torjn- 
stance 390 is an upper bound for d-d. 
|o D y Tarski's [1951] decidability theo- 
rem all spatial stereohedra can be ef- 
ecUvely" "determined. In spi e of this 

r te rssstSntJtS 0 *- 

SSohtdra mentioned in Delaunay 
feldoes not yield a practica algo- 
rithm Among others K°™^veSi 
Koch [1973], a^ /ngel [1981] discovered 

spatial plesiohedra ™^>*\™ cts a 
facets, respectively. tigure iu u * 
plesiohedro, '^^J^Z 

Srobably as interesting as the numbers 
proDaDiy * carefully choosing a 

Lucted via computer [Engel 19811 *o 
excellent sources on til mgs re» 
Grunbaum and Shephard 1^1987^ 

Tt is worthwhile to menxum 
Vo onoi diagrams for regularly ptac ed 
J™ a i so find other applications in 
sites also u numerical 
mathematics^ They apply 

integration [Babenko 19//), to v 
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nvering problems for congruent 

rp^rs 1964], and to statistical 
s [Rogers [Besag 

gations ot lattice sy=" ei nnT ie 
lollison 1977; Conway and Sloane 



■regularly Placed Sites 
now turn our attention to Voronoi 
Z arSng from sets of irregularly 
Tftes A tiling of R d by polyhedra 
sites, a pd the tiling 

d a ceM complex in R it the i s 
7 to face*, that is, each facet of a 
Ion is also a facet of some other „ 

dron in that tiling. 2 ^^nSSf ^ 
by plesiohedra is a cell complex. 

SoSoi diagram 

g f 6 4 'regtns ^ no longer 

tcfltfrise to various cues- 
M f Svfted by the problem offi£ 

5 nsest sphere packW *°Job_ 
posed the following extremal prob 
low big is the smallest possible 
Region with t facets and deftned 
s with minimum distance two? Es- 
f tv,;* and related quantities 
Litd by Mir [1988a, 1988W for 

TwoTk'on general Voronoi dia- 
ls concerned with their combma o- 
•operties and, in parbcu^ with 
ize, that is, their numbers of faces 
3U s dimensions as a function of the 
r n of sites considered. The S iz * o 
noi diagram is an important quan 
«e it relates the amoun t o space 
to store this structure to the m 
ze Answering Crum s prob em, 
ey and Vranch [1977] exhibit an 
ir ily large set of Voronoi polyhedia 
ach pair of which shares a facet 
cample they give shows that the 
a Voronoi diagram in R is lA» ), 
fas also observed by Preparata 
As already mentioned, a Voronoi 
m in the plane has size OW ™xe 
Z es and vertices form a planar 




Voronoi diagrams are related to convex 



araph The size of higher dimensional 
^ P ' t hn* only recently been ana- 
diagrams ta »iy m and Brown 

lvzed completely, ruee ,. in 
59801 showed that Voronoi ^ ^ 
R « are equivalent to certa^ 

^^nStS^rSults on the size 
sense. Hence, Know particular, 

°[ ^Jju^aSTo^ bound the- 

^ TfB Sd 1983]. Exploiting this 
orems L#ronabu;u x ,1 nu m- 

rela^iP.^^^^ 

bers of individual laces 01 . g idel 

Voronoi ^^^^19821 Fig- 

[198 fl m^JEKi original trans- 
ure 11 ■ » lus ^ a ^ the Voronoi diagram 
form that relates the v 
in R 2 to a convex hull m 
stereographical projection. 



a tiling if facet to facet if and only if it is 
£ for faces of any dimensions less than d 
and Ryskov 19871. 



1.2.3 Generalized Voronoi Diagrams 

~i„ r. n + pverv convex polyhe- 
Interestingly, not every ^ ^ & 

dral surface m « mQre gen . 

Voronoi diagram m ft _ corresponde nce 
era l class, a ^^^[1982] and 
can be established J*** 1 showed that 
Aurenhammer [1987 a] sho^ ^ ^ 

Car^olXn^ 
surfaces.* Thta f nerdaa^m a» g 

partiCUl " r nrreSaces the euclidean dis- 
sites p and replaces , 
tance 8(x, p) between a point x ana p y 

" mn „u, railed Schlegel 

No 3, September 1991 
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Figure 12- Apolloniusm 



,2,- v) - w . Intuitively speaking, w ex- 
res^the ~*f^ 

% Center p and radros if ^ > £ 
As an important property, the sepa 
defined by two weighted sites m 
stil l a straight me: th pow« lu*^ 
circles. Figure 4 diagra ms 

diagram in the P l ^ e Diric hlet [18501- 
already appeared m Dirich ' ' d 

therefore, often ^^.^coatinued 

fflirc Tfn est especially in 
to be objects oi ' 1964 i niumin- 

S ?T Kne" ^artWBll. and 
rieomeS of numbers [Gruber and 

Lekkerkerker J^. ^ yes 

S a multipUcatwe constant ^ yield 

Apollomus mo ^^ d ^ sbrunne r [1984]. 
Aurenhammer a? d J d *^"£ the plane 
The separator of two sites mt P 
describes then • ApcUon us circle^ 
Figure 12 it can w partition 

For many additional types, see 
Bolker [1986] hg concept 

At this point, let us roomy 
of Voronoi diagram m another way 
order-fe Voronoi diagram of n sites is 

same ft closest bi complex 
3LS o^n'ed. In the case of fc « 



odel [Aurenhammer and Edelsbrunner 



1984]. 



V- 

L. 
J" 




D1 , QY . or der-2 Voronoi diagram 
Figure 13. Planar oraei a 

[Shamos 1978]. 



i the structure is often referred to 
Lthe furthest site Voronoi diagram smce 
now the region of a site . con 
in space furthest fronr ^igvje 1 ites in 
an order-2 Voronoi diagram for »s 

the Pl r te^tfml aTegU and that a 
pairs of sites detme a 8 edge . 

Voronoi 

-e found m , Mxles [1 JTMgd 
and Hoey [19751 ^ bounds on the 

m P rovl ^ "X rs Tregions, edges, and 
ma ximum numbers 0 Voronoi 

vertices ot a pi<« num - 
diagram. Asymptot a lly t hese.^ 

berS 516 v ( ( belw 'n Voronoi dia- 
relationship between 

crrams weighted aiagrdu^, 
Sams other geometric objects was 
SS out by Edelsbrunner^^ jL 
[1986], Aurenhammer [1987 aj, 
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.hammer and Imai [1988]. Within 
Srt the concepts of power dia- 
zThyperplane arrangement^ 
ra l role Section 3.1 is devoted to a 
3 description of this material. In 
Sax various results on the size and 
"computational complexity of con- 
ing such diagrams are obtained 
e? still little is known on the size 
oroeVk Voronoi diagram in higher 
Ins' consult [Edelsbrunner and 
1986] Only for the furthest-site 
i exact upper bounds are avai - 
Seidel 1987]. The number of re- 
of an order-fe Voronoi (or power) 
* is closely related to the num- 
k-sets of a finite point set. Deter 
g these numbers belongs to the 
open problems in combinatorial 

eThe usefulness of known types of 
oi diagrams had been realized, fur- 
generalizations were attempted, 
ation mostly came from apphca- 
in computational geometry. In or- 
7 meet practical needs, general 
\ e flnJons like the L p -metrics 
l98 0a] or convex metrics [Chew 
rysdale 1985] were taken to . detine 
onoi diagram. Also, the shape of 
,as varied while using a cannomcal 
;ion of the euclidean ^f^ e ^{ 
•jne segments [Kirkpatrick 1979J, 
igments [Yap 1987], and disks or 
d convex sites [Leven and Sham 
19871 have been considered. It trie 
nts form a closed curve, then • dia- 
is just the medial axis of a contour 
ire Figure 9). Disks can be viewed 
ters weighted additively by their 
thus giving rise to the Johnson- 
model. This list is not meant to be 
stive; we shall be concerned with 
matical and computational proper- 
,f varous generalized Voronoi 
ms in later sections, 
us mention one more way of gener- 
on: changing the underlying space. 
h and Im Hof [1979] investigated 
ehavior of Voronoi regions m 
tnn manifolds. Brown [1980], 
nger [1982], and Yap [1987] ob- 
. that Voronoi diagrams on the 



snhere and on the torus, respectively, are 
£ related to their equivalents in the 
euclidean space of the same dimension. 
Diagrams on three-dimensional po y- 
r j i onvfarp^ and on the tnree 
hedral surtaces i._ as ,*. P d bv 

dimensional cone have been treated _ Dy 
Xunt [1985] and by Dehne and Klein 
[1987], respectively. 



1.2.4 Delaunay Triangulations 
Wand in hand with the investigation of 
V^no? diagrams goes the ^sUgation 
nf related constructs. Among them, me 
telSunay triangulation 
nent It contains a (straight-line) eage 
connecting two sites in the P^ne if and 

ZT^y tr^iSm for sites that 
trm a lattice and was extended by 
Eurtay [1934] to irregularly placed 
Stes by means of the empty-circle me- 
thod- Consider all triangles formed by 
he sit- such that the 
Pac h triangle is empty of other sites^ 
The set of edges of these triangles 
Jwes the Dehfunay triangulation of 

th The planar Voronoi diagram and the 
DetaunSr triangulation are duals - a 
orarmtheoretical sense. Voronoi vertices 
Respond to Delaunay triangles 
vX£ regions correspond to > sites and 
Pdffes of both types correspond by oetini 
h ^ Fr om Figure 14 it can be seen that 
SunaTed - (solid) are orthogonal to 
their corresponding Vorono edges 
dasned)-but do not necessarily inter _ 
iect them-and that the boundary of 
the convex hull of the sites consists of 
Delaunay edges. The duality immed, 
ately implies upper bounds of Zn- ,6 i and 
of L - 4 on the number of Delaunay 
edges and triangles, respectively -(com- 
D are Introduction). The Delaunay tri- 
angulation and its duality to Voronoi 
dSams generalize to higher djnjen- 
sions in an obvious way. By results 
by Dewdney and Vranch [1977], ^he 
Delaunay triangulation m R may 
already be the complete graph on 
n sites, thus having (J) edges. For a 
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Figure 14. Voronoi diagram and Delaunay trian- 
gulation are duals. 



catalog of properties of higher dimer, 
sional Delaunay ^angulation ansmg 
from site lattices see Gruber 
Lekkerkerker [1988]. are 
Several interesting properties are 
known for the Delaunay triangulation 
fn the plane. It was first observed by 
Sbson P [1977] that this 
iq locally equiangular, lhis property 
Lis when, for any two triangles whose 
union is a convex quadrilateral the 
nlacement of their common edge by tne 
& diagonal does not increase the 
minimum of the six interior angles con 
^ Actually the Delaunay tnangu- 
Sn is tt oX one with this property 
S particularly shows its unWj 
For sites being in general Position (the 
Delaunay triangulation may contain 
jjeiaunay & cocircu- 

more-sided taces n ioui o» 
iar) Edelsbrunner [1987] showed that lo- 
cal ^angularity is equivalent to global 

SeTSr SUd if only the 
Sest angle of each triangle is consid- 
"eSawTon 1977]. Note that the Dela^ 
nav triangulation thus maximizes the 
Simum angle over all triangulations 
S a given set of sites. On the other hand, 
° a f simple example shows that the j maxi- 
mum angle is not nummaed. I ; is inter 
esting to note that, by the empty circle 
niooerty any triangulation without 
Kse angles" must be Delaunay. Tri- 



angulations without "extreme" an- 
gles are desirable in finite element 
and interpolation methods. 

Lawson [1972] gave a counterexample 
to Shamos and Hoey's [1975] conjecture 
that the Delaunay triangulation has 
minimum total edge length. It does not 
even approximate the shortest tnangula- 
t 10n [Manacher and Zobrist 1979], and m 
fact it may be as long as any triangula- 
tion [Kirkpatrick 1980]. It is, however, 
close to optimal on the average [Lmgas 
1986a]. A different criterion of op > timal 
ity is mentioned in McLain [1976] For 
each triangle, all its pomts should be at 
last as close to one of its defining sites 
as to any other site. This property is not 
shared by the Delaunay triangulation, as 

is claimed there. 

As an important fact, the Delaunay 
triangulation is a supergraph of several 
well known and widely used graphs 
Tanned by a set of sites in the plane. 
Among them are the minimum spannmg 
7™ (or Prim shortest connection net- 
work Produced by Kruskal [1956] and 
Z 1 n Q*S71 the Gabriel graph mtro- 
lucTd by Gabriel and Sokal [1969], and 
tC Jative neighborhood 
rluced by Toussaint [19801 becuon ^ 
gWes definitions of these graphs and 

pk Sed to find certair ; = 

nTcTff co^xTufl^an be rep- 
uted as a weighted mass eerier of th 
■^a o^iarpnt in the triangulation, see 
Sfi^lMW who mentions applications 
£ surface smoothing. Aurenhammer 
1988b] generalized this result to power 
diagram! and their order-fe modifica- 
? ^Tr, reticular, Gale transforms of 
the sites "may be derived in this way 
Surenhamme" 1990b] Gale tranjorm 
are a versatile tool in the investigation ot 
hfgh-dimensional convex polyhedra 

^S'greSd Delaunay triangul.- 
tion obtained from power diagrams gains 
in importance from the follow ng recog 
nit on problem: Given some cell complex, 
? a nTt P be interpreted as a Voronoi di- 
gram? A cell complex in R may 
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, as a power diagram if and only if 
complex admits a certamdual 
m- railed the reciprocal figure oy 
U n r?864] Crapo 11979], Whiteley 
6 and Sh and'Bolker [1986]. The 
W figure is completely character- 
properties of the Delaunay 
Nation for power diagrams Using 
•ocal figures as a criterion 
rmmer S [1987b] showed that sev- 
Sell-known types of cell com 
are power diagrams. Moreover, 
oi di&ram. can be -cognized 
•heir defining sites can be re 
in an efficient way using recip- 
Ggurea [Aurenhammer 1987* For 
p survey on reciprocal figures 
n et 1 [19881. Essentially distinct 
a for cell complexes to be classical 
be multiplicatively or additively 
t b e e d , Voronoi 

hv Ash and Bolker [1985, isboj. 
fi y 967] first suggested that gener- 
haped sites could be -conf uctf 
the shape of the regions , Key de 
Palabi and Hartnett [1968J elaoo 
°on tnis question in some detail 
nition problems of this kind find 
ations, aside from geography ^and 
m ics, in statics and in computer 

-'list of geometric and combina- 
' pSSertuS of Delaunay triangu- 
s is not complete. We shall 
farious others while discussing 
ations in computational geometry. 

Stochastic Properties 



■ we have surveyed a good deal of 
•ch on geometric and combinatorial 
s of Voronoi diagrams. Another 
fcant stream of investigations con- 
the determination of statistical data 
the diagrams obtained l from ran- 
distributions of sites. One _ of the 
'st and strongest motiva ions for 
ing stochastic properties of Voronoi 
Js stems from their practical rele- 
to physical and chemical processes, 
Sly in metallurgy and crystallog- 
iJohnson and Mehl 1939]. Accord^ 
most efforts concentrated on si es 
butedin R 2 and in K . In his valu- 



regi0 ns <**^#ZT»£k as well 
mg tWverage ^number of vertices, edges, 

ances for * e ™ U found some variances 
regions He also io section 
associated with plane or tal 
through region. Ba ed on expe^ ^ 

results, Kiang [19bW gave Voronoi 

regions m R , along w 3 
its generalization to K -d t ^ 

con ecture is, however uu, 
ing to results m ^^ 1 ^ aB in 

fdotfet 0Served ( volume over mean 

volume is scaled. , Wnuen study of 
Miles [1970] nia d e a trough . ^ 

Voronoi of 
homogenous site process^ P imeteT> 
the edge number ed ge lengt 

and area of the ™*° n De l a unay 
Further results ^cern ^ 
triaBgulation a^^J^^ ^ta 

diagrams. R f eva p ?7n 97 8] His ob- 
can be found m Cram W 
served frequencies of edges .per po yg 
are of particular mtere ^™* Mb . 
retical results are [1989 ] 

^ V STt ^hTrs^aTSformly 

mention only ^jgyj^coiu^ 
the^a^^^ 

Lng produced in ^^^abmS 

than being ^ aw ^^ y s ^uch results in 
distribution. We shall see sucn re 

Section 1.3. Drec lude mentioning 

Space constraints pre^ua 

all related work. To name a 1 
we refer to Newman et aL [1983] <n 

ber ° f rwtire td °Sv£r [1984] 

mensions), weaire 
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Figure 



1 5. Empirical density function 



of the volume of regions [Kiang 1966]. 



Swa°y n anT'sioane [1982] (site 

iaUices; applications to spread 0 f 

epidemics and to coding theory are 

cHven) and Aurenhammer et ai. 

fl991] (probabilistic distance ftinctams). 

The reader interested in geometrical 

u v/iit„ in eeneral and stochastic 
probability m geneiai <* . , 

properties of cell complexes in P^** 
mav consult the expository papers by 
Soran [1966, 1969], Little [1974], and 
Baddeley [1977]. 

1.3 Computer Scientist's Viewpoint 

We have documented the remarkable role 
that Voronoi diagrams play m the math- 
ematical and applied natural sciences. 
Y^t for a long time their practical useful- 
ness suffered from the absence of reason- 
ab?y shnple and efficient methods for 
their computation. This section reviews 
me hods for the computer construction 

Algorithmic applications of Voronoi 
diagrams and of related structures are 
discussed in Section 2. 



1.3.1 Early Algorithms 

Since Voronoi diagrams have been used 
for decades by natural scientists, many 



intuitive construction rules were pro- 
posed. The earliest diagrams were drawn 
w°th pencil and ruler; see Horton [1917] 
Ir Kopec [1963] who mentions problems 
of ambiguity if many sites he on a com- 
mon cide. Probably the most obvious 
rploach is to delineate the diagram m 
?he plane region by region, by singling 
out those separators that contributejo 
edses of the current region. A provi 
sional and admittedly inefficient version 
oHuch an algorithm was described by 
burger fl973], Ot^e^X 

rithms build up ^d^K/^ in- 
vertex [Brassel and Reif 979] or by m 
cremental insertion of sites (i e of 
reeions) [Green and Sibson 19 Hi. 

ISure 16 illustrates the insertion of a 
siteT^t solves two tasks. First we 
need to find the current region in which 
Halls. Let q be the site 
Vegion; the separator of p and q then 
will contribute an edge, e, to P s regwn 
Second, the boundary of p's region (bold) 
Ts created edge by edge, starting , with ^ 
During this process, the parts of ^the old 
diagram closest to p (dashed) are 
£££d and deleted. These , par* . aj 
verified by lying on ps siae oi 
he separators of p and its new neig^ 
bors If appropriately implemented, the 
second task requires time proportional to 
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Figure 1 6. Inserting a Voronoi region. 



calculating the dimensional Delaunay 

triangulation but give no concrete 

complexity analysis. 
OtUv recently have optimal or near 
Unly receru j constructing 

optimal algorithms 101 <- 

vUnoi diagr- ^-^^ 

beln fully developed and exploited for 
computational geometry purposes. 



nunber of edges deleted. This num- 
a Od) in the worst case (i denotes 
umber of sites inserted so far) since 
ave to delete a planar graph with at 
i regions. The first task requires 
time in the worst case, but a simple 
urtic will achieve 0(71) expected 
Experimental results showed that 
method is quite efficient and thus is 
cable to rather large sets of sites. As 
the other algorithms mentioned 
e, however, the worst-case behavior 
n 2 ) for n given sites, 
veral intuitive methods are also 
able for the Voronoi diagram m K ■ 
tmpute the induced eeU comp£ 
by facet [Brostow et al 1978], ver 
w vertex [Finney 1979], or via its 
the Delaunay triangulation 
; m ura et al. 1983]. Potentially, they 
ire 0(n 4 ) time in the worst case 
,ugh they may perform quite well tor 
.us distributions of sites, tor d 2 * 
nsions, insertion strategies working 
tly [Bowyer 1981] or based on the 
unay triangulation [Watson 1981] 
been used to construct a suitable 
rinatorial representation oi the 
noi diagram. Their time complexity 
n^ d ) Tand 0(n^ d ), respectively, 
ided the sites are "well distributed. 

should be contrasted, however, 
■ist the maximal size of a Voronoi 



7.3.2 Speeding Up Insertion 



The process of building a Voronoi dia- 
gram in the plane by incremental mser 
STof sites Sands out by i«^ d 
Originally having an O(n^n) fxpeaea 

Ly be polished up to X- 
eral distributions of the sites, oy 
during suitable ordenngs of the sites, trie 
expected time for finding the region the 
expectea i integrating the 

next site tails m anu iu & [Ohva 
npw region can be lowered to U(L) l^ny* 
Z II 1984a 1984b]. 4 This gives an ex- 
pected ruS'me of 0(n). Clearly thu ms 

the best we can b^£^£t3 
of the diagram is a trivial lower u 
for the time needed to compute i 

An alternative approach to speeding 

up insertion is ^dlhor 

nsSing the h sites in random order yields 
a ^ OU log «)-time performance with 
high p obability. This complexity - inde- 
^nt of the Jjtfb-hoa rf OjngJ 

extends to the class called abstract 



^iTn this context that the average number of 
edges of a region is less than six. 
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Voronoi diagrams by Klein £989] with, 
out increase of runtime This general 
concept includes power diagrams and di- 
acn-ams defined by line segments or by 
Srics Guibas et al. [1990] propose 
an even more practical version of ran 
donSed incremental construction of the 
classical Voronoi diagram. 

As with other geometrical algorithms 
the problem of numerical errors arises m 
the construction of Voronoi diagrams. For 
example, sites nearly lying on a common 
drcle define vertices that tend to ap- 
proach arbitrarily close. A PJJ^^ 
making the strategy of insertion robust 
Tgalns! numerical errors is outlined m 
Sueihara and Iri [1988]. . 

The process of inserting Voronoi re- 
signs extends nicely to R ■ Regions are 
fonvex polyhedra that can ^ connoted 
facet bv facet by intersecting existing 
Sons with separators that are planes 

case. A Lgion caiinct 
than n - 1 facets (one for each dnterent 
sS and thus by Euler's relation, has 
0(1) edges and vertices. Therefore one 
needs 0( n + t) time per region if t faceto, 
edges or vertices are deleted during 
its insertion; see Aurenhammer and 
fdelsbrunner' [1984] who treat insertions 
in a more general cell complex m R ■ 
Since'ach Component deleted ^ , be 
constructed first, an 0(n )-time algo 
r Sim results. This is worst-case optima 
since a Voronoi diagram m R may have 
a size of 6(n )• 




1.3.3 Divide and Conquer 
A widely used method to design fast aL 
gorithms is divide and conquer. ^Shamos 
and Hoey [1975] observed that this 
method applies well to problem. jm ^cjm- 
putational geometry and in Particular^ 
So the construction of a Voronoi digram 
in the plane: The given set of n sites is 
dVided into two subsets by a vertical 
£e The diagrams for these subsets are 
computed recursively and are then 
«Xged" in the conquer step to form the 
toUl diagram. Figure 17 shows the ^chain 
of edges (bold) to be constructed during 
1 process of merging the diagram for 



Figure 17. Merging two Voronoi diagrams 
[Shamos 1978]. 



left sites (dashed) with the diagram for 
right sites (dotted). 

Section 3.2 is devoted to a detailed 
description of the method and its gener- 
ations. In particular, two diagrams 
can be merged in O(n) time, which im- 
phes an 0(nlog n) time algorithm if the 
d ivTde step is carried out m a balanced 
wly This is worst-case optimal since any 
algorithm that constructs an explic t de- 
scription of a Voronoi diagram may be 
used to sort: Interpret an input sequence 
of n real numbers as a set of sites on the 
x axis and construct their Voronoi dia- 
X gram Regions of sites for 

Ambers will ^ZSin O^ 
sequence can be obtained m v g 
by scanning through the regons 

Although theoretically fast ana eie 
cant the divide-and-conquer approach 
gfceSSn disadvantage, ^enunta- 
tion details are somewhat j^P 1 ^' 
and numerical errors are likely by 




• Voronoi diagram; ' tSVS^ 
sets of sites. Seidel 19851 P"™" . . x . coor di- 
sites in the plane with re to * e . ^ 
nates takes 0(nlog n) ^^J™ eve 
their Voronoi diagram is part °f th e inp 
unclear whether sorting helps tor vi» 

«r^r Fortune [1988], however, pointed out that 
^ ^Jfthe sites in two different directions low- 
Sr£*3£££l time for constructing the 
diagram in the /^-metric to O( n loglog n). 



Fig 
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ruction: Dividing the plane into 
, w slabs forces the vertices defined 
ps within a slab to approach mfm- 
Joreover, the expected behavior is 
n( n log n). An 0(n)-expected-time 
ithm that applies divide and con- 
to a certain subset of sites was 
ie d by Bentley et al. [1980]. 
ide and conquer is the basis of a 
class of algorithms for computing 
alized Voronoi diagrams m the 
Let us briefly list some of them m 
to give credit to the authors who 
lered them first. Shamos and Hoey 
] reported that their algorithm also 
>s to the furthest site Voronoi dia- 
'in the plane. Hwang [1979], Lee 
Vong [1980], and Lee [1980] consid- 
point sites under the L^-metrics, 
■as even more general distance 
ons were treated by Widmeier et 
)87] Chew and Drysdale [1985], and 
and Wood [1988]. Imai et al. [1985] 
•d sites under the Laguerre dis- 
that is, their power diagram, 
rning the divide-and-conquer con- 
ion of Voronoi diagrams for sites 
general than points, we refer to 
■atrick [1979] (line segments), Lee 
Drysdale [1981] (line segments or 
tersecting discs), Sharir [1985 
rary disks), Leven and Sharir [1987] 
lT convex bodies), Yap [1987] (curve 
mts), and Lee [1982b] (medial axis 
imple polygon). It should be noted 
most of the algorithms just men- 
1 are worst-case optimal. 

Higher Dimensional Embedding 



forming geometrical problems into 
easily understood and solved ones 
an important role in computational 
>,try (e.g., see [Edelsbrunner 1987]). 
i [1979] first perceived the possibil- 
transforming Voronoi diagrams in 
;o convex hulls in R 3 : The sites are 
;d, via stereographical projection, 
joints lying on a sphere (compare 
e 11). This takes O(n) time, 
onvex hull of the resulting three- 
■.sional point set is computed using 
lgorithm by Preparata and Hong 



U977] that requires 0{n log n) time. The 
convex hull is dual to the Voronoi dia- 
gram of the given sites in a geometrical 
sense. Hence O(n) additional time suf- 
fices for deriving the latter. This elegant 
approach matches the optimal 0( n log n) 
time bound, although it suffers from the 
problem of processing three-dimensional 

° b An important feature of the embedding 
method is its easy generalization to 
higher dimensions. Methods for deter- 
mining higher dimensional convex hulls 
are well established [Seidel 1981 1986], 
hence efficient worst-case algorithms for 
computing the d-dimensional Voronoi di- 
agram become available. Their runtime 
however, increases exponentially with d 
according to the maximal size of a dia- 
crram This should be contrasted with a 
recent result by Dwyer [1989], showing 
that O(n) expected time suffices for com- 
puting the Voronoi diagram in constant 
dimensions d for uniformly distributed 

sites , 

Brown's idea was further developed by 
Edelsbrunner et al. [1986] for construct- 
ing the order-* Voronoi diagram family 
and by Aurenhammer [1987a] for con- 
structing the power diagram and its 
order-fe family in general dimensions 
See also Edelsbrunner [1986], Aggarwa 
et al [1989a], and Aurenhammer [1990aJ 
where order-fe diagrams are obtained m 
different ways from their higher dimen- 
sional embeddings. We refrain from any 
details here and refer to Section 3.1 for a 
comprehensive discussion of this mate- 
rial For different approaches to the com- 
putation of order-fe diagrams in the plane 
consult Lee [1982a] and Chazelle and 
Edelsbrunner [1987]. Clarkson [1987] 
speeds up a combined version of these 
approaches using random sampling 
achieving 0(ft« 1+e ) expected time and 



• In fact, since Preparata and Hong construct con^ 
vex hulls by divide and conquer thi . « ijurt 
three-dimensional translation of a divide and 
conquer construction. By Brown's result, simple 
convex hull algorithms [Clarkson and Shor s 1988] 
Lad to simple Voronoi diagram algorithms, 
however. 
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space (for any e > 0) that is nearly 

0P Em?edding in higher ^rnensions yields 
algorithms for various different types 
of diagrams. See Aurenhammer and 
Edelsbmnner [1984] for the mul toph- 
catively weighted Voronoi digram. 
Aggarwal et al. [1989a] for the medial 
aSs of a convex polygon, Edelsbrunner 
et al [1989] for cluster Voronoi dia- 
grams, and Aurenhammer and Stock 
[1988] for the peepers Voronoi 
diagram. 




1.3.5 Plane-Sweep Construction 
Another powerful technique in computa- 
tional geometry is the pl^-sweep tech^ 
nique [Preparata and Shamos 1985] _In 
contrast to the embedding method tins 
technique decreases the dimension of the 
nroblem Intuitively speaking, the static 
'problem of computing a Voronoi diagram 
in the plane is reduced to the dynamic 
"otlemof maintaining the cross section 
of the diagram with a straight line. The 
algorithm simulates sweeping a line 
across the plane from below. At any point 
^time, th'e portion of the dia gi ;am ^w 
the sweep line is complete (Figure ioj. 
Fortune U985, 1987] first observed ^tha 
updates on the sweep line can be imple 
ranted to cost 0(log n) time if a certain 
Sinuous deformation of the diagram .is 
treated. From this deformation the orig 
inal diagram can be constructed m U(n) 
time The method, its application and its 
modification for generalized Voronoi 
diagrams are described in detail in 

The^lane-sweep approach to Voronoi 
diagrams combines simplicity and effi- 
ciency It achieves the optimal 0{n log n) 

me and 0(») space bounds foi -the clas- 
sical type, the additively weighted type 
and the diagram for line segments 
fortune 1987]. The method has been 
successfully applied by Seidel [1988 
t0 Voronoi diagrams with line segment 
obstacles joining sites (and their con 
strained Delaunay triangulations) ^ This 
type had been previously attacked by 

means of rather involved divide-and- 



Figure 18. Maintaining the intersection with a 
line. 



conauer algorithms [Lee and Lin 1986; 
Sngas 1986b; Chew 1989a; and Wang 
and Schubert 1987]. As was reported by 
Rosenberger [1988], plane sweep also 
performs well for planar order-fe dia- 
Jams even when the sites are weighted 
additively. 

L3.6 Delaunay Triangulation Algorithms 
Since the Voronoi diagram and the 
gelaunay triangulation are . duals the 
combinatorial structure of either struc 
t.ure is completely determined from its 
dual Consequently, the Delaunay trian- 
gulation of n sites in the plane can be 
Stained in O(n) time after an 
0(nloga)-time precomputation ot tne 
Voronoi diagram. Several practical ap- 
pliSns, however, solely 
natorial Properties of he ^ Delaunay 

£ ST — fstTe th! 
coordinates of Voronoi vertices. 

Lawson [1972] proposed an alg^thm 
for constructing locally equiangular tn 

angulations by local ^Pf™^ J f 
ing with an arbitrary ^angulation ot 
h! sites edges are "flipped" according 
o th equianVarity criterion (JWed » 

o T , • i 9 £\ until no more sucn ex 
Section 1.Z.4J untu uu nQ771 c UC . 
•i rroe Qr p vpauired. Sibson 115* / M suc 
changes are rec * U1 ^" fh - Droce dure 
ceeded in proving that this proce 




mm 
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19. Equiangularity and empty, circle prop- 



illy constructs the Delaunay tn- 
Llation Lawson's original method 
somewhat speeded up by Lee and 
chter [1980] who proceed mcremen- 
and apply the edge flipping proce- 
after each insertion of a new site, 
lin [1976] uses the empty circle 
>rty to construct the Delaunay trian- 
ion by successively adding triangles 
e circumcircles are empty of sites, 
ct Figure 19 illustrates the equiva- 
of the empty circle property and 
2 quiangularity criterion. Although 
r simple to implement, all these al- 
kms have an 0(n 2 )-time behavior m 
,orst case. On the other hand, Maus 
.] reports that a modification ot 
un's algorithm based on preparti- 
ig exhibits an 0{n) expected run- 
time for uniformly distributed sites. 
> and Schachter [1980] gave an 
og n)-time Delaunay triangulation 
ithm; it uses divide and conquer and 
oilar to Shamos and Hoey's primal 
loi algorithm. A modified rniple- 
ation using prepartitioning was 
n by Dwyer [1987] to run m 
oglogn) expected time while re- 
ig the optimal 0(n log n) worst-case 
it generalizes to the L p -metrics for 
I. Drysdale [1990] modifies the 
3-and-conquer approach to compute, 
• n log n) time, Delaunay tnangula- 
' generated by general convex dis- 
functions. A recent algorithm by 
is et al. [1990] inserts (point) sites 
idom order. Using a structure simi- 
. the Delaunay tree [Boissonnat and 
md 1986], earlier versions of the 



triangulation are maintained to facili- 
tate the determination of the triangle 
covering the next site to be inserted. This 
practical algorithm exhibits a random- 
ized running time of 0(»log ») and lis 
similar in spirit to Clarkson and Shor s 
[1988] convex hull algorithm. 

The question arises as to which of the 
planar Delaunay algorithms can b€ ^ ex- 
tended to higher dimensions. We have 
"already seen some early (althoug .widely 
used) algorithms at the beginning of this 
section. It should be further observed that 
the method of higher dimensional em- 
bedding actually first constructs the 
Delaunay triangulation (via a convex 
hull), then derives the Voronoi dia- 
gram via dualization. Recently, Rajan 
Fl9911 proved a property ol a- 
dimensional Delaunay triangulations re- 
lated to equiangularity. It implies an 
incremental algorithm that makes the 
triangulation locally Delaunay after each 
insertion of a new site by applying a 
triangle-flipping procedure similar to 
Lawson's method in the plane^ The algo- 
rithm is worst-case optimal for odd d. 
Whether every d-dimensional triangula- 
tion can be made Delaunay by local 
improvement remains open. 



7.3.7 Storage Representation and Dynamization 

The Voronoi diagram may be viewed as a 
data structure that organizes its defining 
sites in a prescribed manner. A host of 
applications of this data structure m 
computational geometry will be given n 
Section 2. Out of a number of possibili- 
Ses to store the planar Voronoi diagram 
or, equivalents, its dual, the Delaunay 
triangulation, the ^ad-edge structure 
proposed by Guibas and Stolfi [1985] is 
narticularly practical. Essentially, each 
edge's stored P as a pair of directed edges. 
Each directed edge, in turn stores the 
vertex it originates from and pointers to 
the previous and to the next (directed) 
edge of the region to its left. The quad- 
edge data structure provides a clean 
separation between topological and geo- 
metrical aspects and supports the imple- 
mentation of standard techniques for 

ACM Computing Surveys, Vol. 23, No. 3, September 1991 






366 



Franz Aurenhammer 



computing Voronoi diagrams such as d 
vide and conquer or ^mental inse 
lion. The structure has ^e^^^d 

f rt Wrier dimensions by DODKin ^ 
1^1989] (facet-edge gJ^JJ 
p3^ and by Brisson [iy»yj ^ ceil 
# ) ana These data struc- 

tuple structure m R ) • Th *se a 
tures represent the inciae 

ordering -^—^ See also 

cell complexes th Voronoi 

Like any data structure 
j- „„„ fan be dynamized, tnat it., 

dia ^f „^ for a set of sites that varies 
maintained for a set o ^ 

over time by CQVered by the 

case of site insertion , that build 

^t^Wf™ of ? £f» 

into the Voronoi diagram for n sites in 

tions and deletions « y or onoi-free 

tree records the history 01 diagram 
conquer construction of the dia f in 
rPianre 20) Its leaves hold the sites m 
(figure lu). tt n rh inner node is 

5223 ^h^ r dSam n r comes 
JroTcombining its sons' subdiagram. 
SETidea originates with Overmar [1981] 
who, however, did not give details oi the 
Z hod or its storage requirernen ■ By a 
result by Aggarwal et al 

^Usinfthe Delaunay triangulation a 
plan^ set of sites can be organized ante , a 

Xch a given point falls m efficient 
ovnpoted and worst-case time, me 
Sufay tree reflects the h story 
of an incremental construction of the 



Delaunay triangulation (Figure 21). 
Sny Sngle T destroyed by insert- 
£Sf a new site gets as sons (solid 
^vL+orO the new triangles sharing 
an edge with T . Remaining triangles ad- 

Jr,+ tn T in such an edge get the re- ^ 
Ja >■ ,?Jw triangle as a stepson (dotted ' 
SSSrSS^et al. [1990] further 
Sped this structu- 
results, they showed that the expected 
number of structural changes > of ti» 
Delaunay tree is O(n) if n sites are 
Serted in random order into an (im- 
Slv empty) Delaunay triangulation. 

A different way of dynamizing Voronoi * 
diagrams or Delaunay triangulations is 
to Slow movement of sites (continuous 
updS rather than the r m-arbon* 
deletion (discrete updates). The Voronoi „ 
diat am of continuously moving sites will 
change its shape continuously, although 

^ of thTdt^ 

^ ^ZTte S ener f e cases, 
Sch combinatorial changes always are 
?ocal Sure 22). An edge between two 
„« rwe 2 and 4) collapses to a ver- 

final environment, this procebb 
dual env fe diagonals of the 

angles This sugge moving 
the DeZaunay history °\ a set 0 d u t of 
• tK the chronologically ordered u» 

edg in thek DelaUnay frTthe 

Hon Tokuyama [1988] considered the 

s jLoWdiy ^ 

m ^it-ps respectively. i" BU , 
m, srteb, Mnrri) and can be 

his tory has length 0™)^^ in 

computed ana y v re t r ieve 
0(nm log nm) time m or der ^ to ret 
the Delaunay triangulation and^ 
Voronoi diagram) m 0(n + g JO 
anv eiven moment. For n .sites i 
any give : ndivic iual directions ana 
in fixed but lnaivmu vplocitv the 
with constant individual velocity , 
Delaunay history has length CA J 
This result by Imai et a \ [198 ^ t ^ tive 
be contrasted with the intui 
argument that (j) cocircularities may 
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Voronoi tree [Gowda et al. 1983]. 



r among n sites. Guibas et al. [1991] 
red that the bound remains nearly 
- if the sites are moving contmu- 
) along rather general trajectories^ 
also Aonuma et al. [1990] for related 
Its. 

: Computing Voronoi Diagrams in Parallel 

.llelizing algorithms in computa- 
il geometry usually is a complicated 



task since many of the techniques used 
(incremental insertion or plane sweep, 
orTns'ance) seem inherently sequela . 
Since it is one of the fundamenta struc- 
bmce ii is o Voronoi diagram 

tures m this area, tnc >v 
ha? been among the first geometrical 
structures whose 8 construction has been 
parallelized. The underlying model ^ 
computation mostly ^as been the CREW 
PRAM (concurrent-read, e^^^J 
parallel random access machine), where 
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Figure 21 . Delaunay tree. 





Figure 22. Combinatorial change. 



processors share a common memory 
that allows concurrent reads but no 
simultaneous writes to the same cell. 

As has been mentioned, fl(nlog n) is a 
lower time bound for sequentially con- 
structing the Voronoi diagram for n sites 
£ the plane. Therefore, 0(log n) parallel 
time is the best that can be achieved 
when only O(n) processors are to be used. 



In her pioneering work, Chow [1980] 
showed that convex hulls in three space 
can be computed in 0(log 3 n) time usmg 
Oin) processors. By their duality to 
planar Voronoi diagrams the 
carries over. Aggarwal et al. [1988] im- 
proved the Voronoi diagram construction 
to Oflog 2 n) parallel time by applying 
dWide and conquer. Their approach criti- 
cally depends on restricting the domain 
Sat can contain the "merge chain 
(compare Figure 17) and has beer uq- 
proved recently by Cole et a . [1990] 
using a tailor-made data structure. An 
optimal parallel construction of Voronoi 
diagrams is still outstanding, however 

Goodrich et al. [1989] construct the 
Voronoi diagram for line segment sites 
in CKlog 2 ) time using O(n) processors. 
Their approach to constructing this im- 
portant generalization is similar in ^spint 
to Yau's [1987] curve segment algoritnnv 
As a related question, Schwarzkopf 
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r iQ89l addresses the problem of comput- 
er the digitized image of a planar 
Voronoi diagram using certain processor 
n etworks. 

2 . ALGORITHMIC APPLICATIONS 

The precomputation of a Voronoi dia- 
gram is the initial step of various algo- 
rithms in computational geometry. We 
have already mentioned some of them in 
the Introduction. In this part of the 
survey we demonstrate the broad 
scope of algorithmic applications ot 
Voronoi diagrams and of closely related 
structures. 

2.1 Closest-Site Problems 

According to their definition, Voronoi di- 
agrams apply naturally to various prox- 
imity problems. Determining closest sites 
plays a dominant role in this context 
since it appears as a subroutine of many 
geometric algorithms. Applications in- 
clude clustering and contouring and var- 
ious other problems whose relation to 
proximity may not be so obvious. 

2.1.1 Nearest-Neighbor Queries 



Probably the most popular problem in 
this area is the post-office problem. How 
should a fixed set of n sites in the plane 
be preprocessed in order to determine 
quickly the site closest to an arbitrarily 
chosen point (the query point)? The post- 
office problem is motivated in trie 
Introduction as a two-dimensional 
file-searching problem. Another impor- 
tant application stems from a common 
data classification rule that requires as- 
signment of new data points to the same 
class in which their nearest neighbor lies. 
An early solution yielding 0(log n) query 
time but requiring 0(n 2 ) storage was 
provided by Dobkin and Lipton [197bJ. 
Shamos [1975] pointed out that the 
Voronoi diagram of the sites divides the 
Diane into regions of equal answer with 
respect to the post-office problem: By del- 
ation, the region of a site contains all 
joints closer to this site than to all other 



sites This reduces the problem to finding 
the region that contains the query point, 
an approach that became powerful once 
efficient data structures were developed 
for the latter problem: point location m 
planar straight-line graphs Kirkpatrick 

1983; Edahiro et ^J^nX 
Edelsbrunner et al 1986L », 
logarithmic query time and linear 
stoVage are achieved. This ; is asym^ 
totically optimal since it matcnes 
the information-theoretical lower bound 
Shamos' approach in conjunction with 
Kirkpatrick's was worked out in a geo- 
metrically dual setting by Edelsbrunner 
and Maurer [1985]. Usage of a structure 
related to Boissonnat and Teillauds 
U986] Delaunay tree even obviates tne 
need of postprocessing the Voronoi dia- 
gram for point location; Gmbas et al. 
[1990] show that post-office queries are 
supported in 0(log 2 n) expected time by 
this structure. An off-line solution to the 
post-office problem (all query points are 
given in advance rather than being spec- 
ified on line) was presented by Lee and 

Yang [1979]. . . 

In principle, higher dimensional vari- 
ants of the post-office problem may be 
solved via the Voronoi diagram too. 
Their practical relevance is evident irom 
the equivalent file-searching problem for 
multiattribute data mentioned m the In- 
troduction. Although Chazelle [1985 
showed that point-location in a Voronoi 
diagram in R 3 can be performed effi- 
ciently, the approach suffers from a con- 
siderable storage overhead of 9(/i ) that 
comes from the worst-case size of the dia- 
gram In d>3 dimensions, the pro- 
hibitively large size of the diagram may 
be reduced to some extent for that 
purpose [Dewdney 1977]. 

Generalizations of the post-office prob- 
lem for nonpoint sites or for modified 
distance functions may be solved by 
means of an appropriate Voronoi dia- 
gram as well. We particularly mention 
Mount [1986] (geodesic distance on a 
polyhedral surface) and Aurenhammer 
et al [19911 (probabilistic distances in 
the plane) as two generalizations leading 
to interesting types of diagrams. 
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2 1 2 k-Nearest-Neighbor Queries 

An obvious and ^^erf^ 
of the post-office joKemjs to ask ^ ^ 
k closest sites for a q"«g^ fl975] that 
observed by Shamos and Hoey 119/ J 
the order-fe Voronoi diagram iOi 
plane into ^^leTp^S^B 

the diagram ioi . f handles queries m 
data str ^ re tim h e a l n d a ^ 0(fe 2 n) 
0 (log n + « ma ^ 
storage, line a»8 hich 
0(fc(n-fe)) regies for eacn o 
th e k closest sites are stored ^1 
specified in the input an U(n ) p 

daU ^oUaUd l^Lrner't af 
time was obtained by r, u rf 

[1986] via r embedding ^ 

0 rder-fe Voronoi diagr ams io 

into an arrangement of ^pl { di . 

2T?lSSd ^astically'educes the sto. 
%fr^rL t of these Jg^ft 

<r£ J P 1oble P m that 
all sites covered by a given juery 
In facility location for exampl 
mig ht be interested m all sites ^ 
influenced by a new feed ity oi g & 

dius of ^ tra t?H on F we1e undertaken by 
satisfactory solution w ere^ un 

Bentley andMaurer order . fe 
is based on a certain ^uence ot ^ 
Voronoi diagrams of the give 

has been refined by Chaz ^ al et al . 
[1986b] and recently by Affi^** 
riQQm The latter refinement y ieia& 

g£ts» - to, &r s 




Figure 23. All closest pairs and the largest empty 
circle. 



2.1.3 Closest Pairs 



Z »™l£ty £oK*™- 0ne obvious ff 

3 Arises in collision detection; the 
?^o dosest Ss are in the greatest dan- 



ger of f*2^£^& 

example, in clustering o g . tes ^ 
the closest for eacho t 

Strai t tf Ef eC^stances" In the 
in calculatm Lh lotions were settled 
planar case both questio ^ 

in A im .v nT75] b ^ prlcomputing the 

23 r^^s^ Seg- 
ments. Since ^ he regi gh 

region of its closest s ed 
edge, one only needs to vi s 

f nd Ca ! C e : la whoi e separator Yields that 

two sites wnose »p 

edge . This clearly t^ O^) *m 
the diagram - availa ble^ ^ 

K Offices tJ Sain certain parts 
that it sumces sweep 

of the diagram during * J\. ex . 

rather than to compute the .^g^ 

pl ici tly in ordeMo-lv both probl ^ 

optimal time. Curiously, w / 

applications. 
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22 placement and Motion Planning 

v„rr,noi diagrams have the property that 
V °TZnt on a Voronoi edge maximizes 
^ Stance to the closest sites. Al- 
,2h being straightforward, this prop- 
th ° U t essential for the usefulness of 
Voronoi diagrams in placement and 
motion planning problems. 

2 2. 1 Largest Empty Figures 
The largest empty figure problem is the 
UlowX Given n point sites in the 
tZ where should a figure of pre- 
SKd iape be placed so that its area^s 
maximized but no site is covered? If the 
Kreis a circle, the problem can be 
Emulated as finding a new site being 
as toas possible from n existing ones. 
Shamos and Hoey [1975] mentioned the 
reliance of this facility location problem 
operations research and ^tnal e»j 
Serine They also showed how to find 
a Sesfempty circle in 0(n\o g n) time 
Tie Inter of such a .circle is erther a 
vertex of the Voronoi diagram of the sites 
o the intersection of a Voronoi edge and 
the boundary of the convex hull of : the 
sites; see Figure 23 illustrating the oc 
currence of the latter case. Thus atten 
Z I can be restricted to 0(») possible 
placements of the circle. 
H The smaUest enciosmg circle of the 
sites is obtainable in O^ n)*™™ * 
similar manner. Sliaxnos ^l 97 ^ ^"J^ 
that this circle always will be centered at 
a vertex of the furthest site Voronoi dia- 
eram Tins was proved true by Toussamt 
fnd Bhattacharya [1981]. From Figure 
24 it can be seen that the edge ; of this 
diagram form a tree structure _ This re- 
stricts the maximum number of vertices 
and thus of placements of the circle to 
n - 2. 7 On both problems, there exists 




Figure 24. Furthest site Voronoi diagram and 
smallest enclosing circle. 



7 It is easily seen that the regions of the furthest 
site diagram are either unbounded or empty . The 
**Sg furthest site Delaurmy tnangulatwn 
Zfan outer' planar graph (with at most » ^ 
as vertices), which is well known to contain at most 
2n - 3 edges and n - 2 triangles. 



a considerable literature 
and Shamos 1985] In P^f^^ 
problems are solvable m O(n) time by 
Lans of linear programmmg as was 
demonstrated by Megiddo [1983L 

The Voronoi diagram approach gener 
alizes-while still achieving the 
oSff n)-time bound-to other figures 
these are "circles" with .respec 
to the underlying distance function^ For 
example, the largest empty axi^parane 

square can Z 
rliapram m the i^-mew."*- ^ 
S^Se [1985] find largest empty 
nomothetic copies of 'genera Uy .shaped 
convex figures via diagrams that are 
induced by a rather general class ot 
distance functions. The 
the largest empty axis-parallel rectangle 
turns out to be somewhat more compli- 
caS the aspect ratio of the rectangle, 
Sd thus the distance function to be used, 
L n* known in advance. This problem 
has an obvious application. If we are 
^Tven a piece of fabric or sh eet ^etal and 
the sites are flaws, what is the largest 
area rectangular piece that can be sal 
vaeed? Chazelle et al. [1986a] gave an 
0(n log 3 ^)-time solution. They use an 
Tntereling type of Voronoi diagij^jt 
can be intepreted as the power ^ diagram 
of a certain set of circles [ Edelsbr , un ^ 
and Seidel 1986; Aurenhammer and Irnai 
19881 The running time above was 
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Sun [19871 « a di agram for line 

Pr ^^^rbin exploited by 
segment sites n* largest 

Ch TXp^ ^3 -vex polygon such 
similar copy 01 d As a re . 

that no segment is overly 
lated question, AonumaetaLUy^iP 

a given ^-^g^iZdLy of an 
far as possible from "ie relies 

arbitrary P^^Ims dySmized with 

nn Voronoi diagrams) 

Aspect to certain site movements. 



I / 



2 2 2 Translations! Motions 
Voronoi diagramsare an aid ir i the plan, 

when moving on *e edges o 

let™. » wjcrxs 

Rowat [19791; several authors p 
this idea father We to 

For /W^S) of obstacles by 
approximate the scena ^ Yap 

polygonal sets, w between 

fl985] plan ^^J^^^ut. In 
two given pomts m this e™ : ^ 

5 *eir .-^yrsS 

Sagram can be computed L > 
time [Kirkpatrick 1970, tor 
Figure 25 depicts^ the g^apn ( ^ 
four line segment obstacles ( > 
serve that each point on an eo-B 

•jiatant from its two closest obsta 
T Therefore for each edge, it is easy 
cles. lnereiore, ^ pipnrance when 

to compute ^ T ^^SKU*» 
movm g on i • J ^ first fmds the 
motion (if it f™ tQ the C enter of the 



Figure 25- Voronoi diagram for line segments. 

the clearance falls short of the disk 
radius. All of that can be done in 

^Comp^e^to Thers, this simple ap- 
prSchhls the additional advantage that 
C disk radius need not be known dur- 

S^XTguve a do for th, 

robot as well. 



2 2 3 Rotational Motion 

Son has been P^^— 

amidSt t ?haf is aUed to rotate. The 
Se ^fon o?a segment (of given length) in 
?Wan ^ determined by the triple p = 
x ^ indicating the coordinates of one 

indent and the ^^^K 
all triples p such that the ^ ^ 
segment L(p) does n o underlying 
obstacle is a subset oi . ce . 
three-dimensional J^^^e a 
Within this subset, we may 
Voronoi diagram as the set of all P 

that L(p) i % e ^ 1 t n l a n tss is with re- 
closest obstacles; ^^J^ betwee n 
spect to the minimum distance d 
segment and obstacle 

O'Dunlamg et al. \ l ™°> £ .,. J ia . 
■aJa » thorough analysis of this aw 

;£ {LS 4r - rinV— 
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♦meted and a motion be planned in 
TnTniol* n) time, thus being a sig- 
°V\nt improvement over earlier tech- 
°5S sSSy and Sharir [1986] speeded 
oiq ^f method-by introducing a related 
UP fl th^o run in 0(t log n) time where 

%lSse^ S£ ^ons of the 

^xlomati ^ chSterization of the 
wndar P roperties of a generally shaped 
SSnsfonal scenaro in order -to £■ 

gfa Voronoi d^.^^SSSf 
I!l«nninB was proposed by Stifter uvovi. 
£ a more detailed introduction and re- 
, nT tbe use of Voronoi diagrams in 
Son' pUnSg, the interested reader 
Tay consult Schwartz and Yap [1986] or 
Alt and Yap [1990]. 

2.2.4 Path Planning 




A special case of planning collision-free 
Ltions, where the robot is ******* 
sLgle moving point, is path P^«Jf- 
Finding euclidean shortest paths m the 
p esence of polygonal («r PO^JraB ^ ob- 
stacles has received particular interest. 
For example, consider the problem of 
determining that location for each Mrf m 
fixed locations within a Polygonal factory 
floor, reachable on the f ortest path^ This 
nmhlpm can be solved m U(n + 
Sog V + m) time by first constructing 
fe Voronoi diagram of n sites under the 
geodesic metric interior to a simple pgy- 
gon with m edges [Aronov 1989L The 
site yielding the shortest path to a par- 
£lar site defines a neighboring region 
in the diagram. This improve .a* earlier 
result by Asano and Asano [1986]. Ke 
J ^questions can be solved by means of 
the furthest site geodesic Voronoi dia 
gram inside a polygon; Aronov et al 
[1988] showed how to construct this 
structure efficiently. 

The following important question is 
stm unsolved. Is it possible to compute 
the geodesic Voronoi diagram for n sites 



Figure 26. Geodesic separator. 

and a collection of disjoint and con- 
vex polygons with a total of m edges in 
0(71 + m)log(n + m) time? The mam 
problem in this respect is that com 
nutine geodesic distances is not a 
oti- ime operation. Figure ,26 shows 
he separator (dashed) of two sites under 
the geodesic distance among two convex 
polygons (bold). It is composed of straight 
Kind hyperbola 

intersections with visibility lines (solid) 
Path-planning problems m three 
dimensional polyhedral spa^veto 
studied, by among other authors Mount 
U986], Mitchell et al. [1987 j, ana 
Baltsan and Sharir [1988]. For example 
finding shortest paths lying entirely on a 
[generally nonconvex) polyhedral surface 
interest ^J^SS 

TtaSby using Voronoi diagra^ 
fined in such spaces. Moreover, since 
finding shortest paths is closely related 
to detuning visibility bet^jg**; 
cles, Voronoi diagrams involving visibil 
Z constraints have been exploited in 
this context. An example 1S the peeper s 
Voronoi diagram used by Baltsan^nd 
Sharir [1988] and by Aronov [198^^ 
point in the plane is assigned to the clos 
est visible site, and visibility is con 
tainea to a segment on a Uneavoidu* 
the convex hull of the sites gee Figure 
36). This structure can attain a size ot 
0(V) and is constructable in On ) time 
[Aurenhammer and Stockl 1988]. 



2.3 Triangulating Sites 

Intuitively speaking, trian t ^ l ^ t ^ e 
sites in the plane means partitioning the 
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plane into triangles whose vertices are 
all and only the sites. Triangulations are 
important because they (as being planar 
graphs) on the one hand have only linear 
size and, on the other hand, establish 
connectivity information among the sites 
that suffices for many applications. 

2.3. 1 Equiangular Triangulations 
An early and major application is the 
interpolation of functions of two vari- 
ables, where function values initially are 
known only at irregularly placed sites. 
Practical evidence stems from finite 
element methods or from processing 
graphical data, for example, in terrain 
modeling. Given a triangu ation of 
the sites, the function value at an arbi- 
trary point can be computed by interpo- 
lation within the triangle containing that 
point. Lawson [1972] and McLain [1976 
have reported that a triangulation is well 
suited to interpolation if its triangles are 
nearly equiangular. The Delaunay tnan- 
gulation is the unique triangulation that 
is optimum in this sense [Sibson 1977L 
Sibson [1980] further suggests the use ot 
local coordinates obtained from sites be- 
ing adjacent in this triangulation m or- 
der to compute smooth surfaces. There 
might, however, occur more-sided poly- 
gons (rather than triangles) if various 
sites are cocircular. Mount and Saalfeld 
[1988] showed how to triangulate such 
polygons in order to retain equiangular- 
ity This process does not affect the U(n 
log n) time needed for constructing the 
Delaunay triangulation. 

Given a set of sites and a particular 
function value (height) at each site any 
triangulation of the sites defines a trian- 
gular surface in space. The "roughness 
of such a surface T may be measured by 
the Sobolev seminorm 

£ |a|(«2 + 02), 

where I A denotes the area of the spa- 
tial triangle A, and a A and ^ denote the 
slopes of the plane containing A. By ex- 
ploiting equiangularity, Rippa [1990J 
showed that the Sobolev seminorm is 





Figure 27. Constrained Delaunay triangulation. 

minimized if the underlying triangula- 
tion of the sites is Delaunay. This re- 
sult is somewhat surprising since the 
Delaunay triangulation- although itself 
clearly being independent of the heigh 
at each site -optimizes a quantity that 
depends on these heights. 

2.3.2 Constrained Triangulations 

In a constrained Delaunay triangulation 
prescribed edges are forced m as part oi 
the triangulation. See Section 3.3 for a 
formal definition. Figure 27 depicts such 
a triangulation with two prescribed edges 
(bold) This structure can be constructed 
in O n log n) time [Seidel 1988] and has 
two major applications. First, it provides 
a more realistic approach to modeling 
terrain surfaces since salient elements 
like mountain ridges or valleys may be 
prescribed. Lee and Lin [1986] showed 
that the equiangularity criterion is still 
fulfilled by constrained Delaunay trian- 
gulations. Since triangles with small an- 
gles produce a poor computer graphics 
display, the best possible visualization is 
achieved. Triangulations inside a simple 
polygon are an important special case 
The boundary edges of the polygon will 
not be edges of the (unconstrained) 
Delaunay triangulation of the polygon 
vertices, in general, and thus have 
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f be prescribed. 8 Chew [1989b] 
l d constrained Delaunay tnangula- 
"*L q to generate triangular meshes 
Sde polygons where all angles are 
tit ween 30 and 120 . 

Second, constrained Delaunay triangu- 
, 2™, are an aid in path planning. In 
^ to plan a path from A to B that 
Sds a Vet of n line segment obstacles 
S thus arbitrarily shaped polygonal 
regions, for instance) one might con- 
Set this triangulation-by taking A 
and B and the segment endpoints as sites 
and the segments as Prescribed 
Sees- and then search the graph formed 
D v triangulation edges. This was ob- 
served by Chew [1986], who also proved 
following interesting result: Between 
any two sites there is a path m this 
^angulation whose length is at most 
VlO times the geodesic distance between 
them provided the L.-metric (instead of 
ne^udidean) is taken to define the tri- 
angulation. This yields an 0( ? log £ 
time algorithm for approximating the 
optimal path, which is a significant 
improvement over exact methods. 

Dobkin et al. [1990] proved ^a factor of 
about 5 for (nonconstramed) euchdean 
Delaunay triangulations and the eu- 
clidean distance. Their result has been 
improved to a factor of about * by Keil 
and Gutwin [19891. Results of this kind 
are rather surprising since they show the 
existence of a sparse graph being almost 
as "good" as the complete graph on n 
sites, independently of n. 

Let us mention another application. 
Given n sites and some prescribed (non- 
crossing) edges between them, a diagonal 
is a new edge that does not cross any 
prescribed one. A popular heuristic for 
triangulations with minimum total edge 
length is the greedy method. This method 



iteratively adds shortest diagonals by 
considering the current (and ^initially 
empty) set of edges as prescribed^ Lmgas 
[1989] reports that a shortest diagonal 
can be found in 0(») time given the 
constrained Delaunay triangulation. This 
gives an 0(n 2 log ;i)-time and 0(n)-s^ce 
Algorithm for constructing such greedy 
Triangulations. Previous methods re- 
quired 0(n 2 ) space. Whether a minimum 
length triangulation for n sites in the 
plane can be computed in polynomial 
time is an important open question. 



2.3.3 3D Triangulations 



8 It is an interesting question how many «w«Jm 
have to be placed on the polygon boundary in order 
^ force it in as a subgraph of the Delaunay tnan- 
gulation. Schwarzkopf (personal communication 
W90) observed that n(»*) sites may be necessary 
for an vertex polygon. A known upper bound I is 
2 0(r.lo gn ) xhe O(n) bound claimed in Boissonnat 
[19881 involves a constant depending on the angles 
of the polygon. 



The problem of "triangulating point 
sites in three-space arises in the decom- 
position of three-dimensional objects for 
example, of computer-generated solid 
models or of objects specified by planar 
Sots sections. Field [1986], Boissonnat 
[1988], and others report on advantages 
of tetrahedral decompositions and pro- 
pose the use of the Delaunay ^angula- 
tion. Although there is so ™ evidence .for 
the quality of thus produced tetrahedra, 
some important properties of the two- 
dimensional Delaunay triangulation do 

^SrtheTumber of tetrahedra may be 
B(n 2 ), 'hence leading to a storage re- 
quirement that is Prohibitive for many 
applications. This fact suggests the use 
of postprocessing. Chazelle et al. [1990] 
showed the following somewhat counter- 
intuitive result: The number of tetrahe- 
dra may always be lowered to about rc V n 
by adding about -Jn new sites. This pro- 
cess will also eliminate elongated (and 
thus undesirable) tetrahedra. In act, 
adding O(n) (well-chosen) sites will al- 
ways f educe the size of the triangu ation 
to O(n), with the additional effect that 
only constantly many edges emanate 
from each site [Bern et al. 1990J. 

Second, the question of finding a 
three-dimensional analog to the eqman- 
gularity property is still unsettled. Rajan 
[1991], however, recently succeeded .in 
proving a related property. The three 
dimensional Delaunay triangulation 
minimizes the largest containment 
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radius of the tetrahedra; that is the .ra- 
dius of the smallest sphere containing a 
Srahedron. This sphere comcides ; w th 
the circumsphere if and only if its .center 
lies within the tetrahedron. Note that 
2 containment radius 
trahedron size more appropriately than 
the circumradius since J^jSy 
eated tetrahedra may have arbitrarily 
farge circumradii. The Delaunay trian- 
2tion is the most compact tnangu- 
£! in this sense. This . property is 
auite remarkable as-unhke the two 
imensional rit«ation-the rfte- 

trah edra ^^S^^S 
triangulating the given si 
vary from linear to quadratic 

There is another property of ^aunay 

fZtZ may be defined with reject to > * 
fixed viewpoint. Generalizing a result by 
De Floriani et al. [1988], Edelsbrunner 
[1989 proved that this relation is acyclic 
for the triangles in a three-dimensional 
Delaunay Angulation, - £*-£S 

™hi£ Sat eliminates hidden surfaces 
ordering the thre.dimensiona 
obiects with respect to the in-front/ 
behind relation and then drawing them 

so-called a-shapes are made of such 
triangles, a-shapes have been ^ujed J 
Edelsbrunner et al. U^l 10 IUU 
shape of point sets. 

2.4 Connectivity Graphs for Sites 

The Delaunay triangulation contains as 
Subgraphs, ^ZT^tX 

sequel. 



straieht-line connection of the sites with 
mSum total edge length. See Figure 8 
for an illustration. Not surprisingly, this 
structure plays an important role m 
transportation problems pattern recog- 
nition and cluster analysis. Construc- 
tion 'methods working on general 

& -d/more recently Yao U975]. 
whose algorithm works in 0(m log log n ) 
time for a graph with m edges. All these 
algorithms use a basic fact: For any par- 
tition of the sites into two subsets, the 
shortest edge (edge of minimum weight) 
Stween the subsets will be present in 
the tree. Applying Yao' s algorithm di- 
rectly to the complete euclidean distance 
Saph for the sites results m an 0(n 2 ) 
running time, however. 

Shamos and Hoey [19751 recognized 

that the edges of a 

tree must be Delaunay edges. This inter- 
esting and important property holds 
See by the fact mentioned above, each 
tree edge is the diameter of a disk empty 
!Hit7s and by definition, each triangu- 
St To edge is a" chord of some disk empty 
of sites From the computational yiew- 
r v/n's algorithm just needs to be 

suffices to derive triangulation 
tree from the ueiaun *J & 
fPreparata and Shamos 1985]. This gives 
an OUlog »)-time algorithm, which i 
Stffl by reduction to sorting n real 

^relationship above extends to more 
general metrics, 

ri 9791 extends it to the L x -metric (.com 
nare FiSe 3). One may be tempted to 

the sites (the dual of the furthest 



i 

i 




ilsiSIiilil 




2.4.1 Spanning Trees 

Given n point sites in the plane, a 
feuclidean) * minimum spanning tree is a 
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*The minimum spanning t-^ould no be rf 

fused with the ^.SSUpoW* 
sites, a concept where the .^""^ge length, 
allowed in order to nun«ni» *e total edg 
The problem of constructing a ™ mmuin 
tree is NP-complete [Garey et al. 1976]. 
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, ni rliaeram). Shamos [1978] showed 

^Tthls coSecture is erroneous 

£ rfpin [1990], however, pointed out 

K two rte. are connected in one of 

S k smallest spanning trees only i 

u % belong to the same subset of sites 

TLinl a region of the order-(fe + D 
defining a r g dg tQ an 

SSSd for finding the * 

^S^ion of the m ini g 
™ spanning tree in three space has 
SSdKta a difficult Problem. Pre- 
mutation of the Delaunay triangula- 
on does not help in general since he 
Imnlete graph may be the output m the 
vTst casf Efficient approximation and 
t d-time algorithms were given 
among other authors, by Vaidya [1988aj 
and Clarkson [1989]. The ^st algorithm 
mnning in subquadratic worst-case time 
Sue to Yao [1982]; it can be speeded up 
o run in 0(»log*) 15 timely using 
to run I" v „ ° n f rarefullv chosen 
Voronoi diagrams ot careiuny 
Groups of sites. By exploitmg-m addi- 
Sn-an interesting relation to a certain 
L P0 i n t problem, Agarwal et al. 
So] were able to bring down the com- 
mexity to 0(n 4/3 log 4/3 n). It still re- 
E open whether an algorithm 
ma ching the best currently known lower 
bouS of ? 0(» log n) can be developed. 

2.4.2 Spanning Cycles 




A graph connecting n sites is called 
Hamillnian if it contains some cycle 
passing through all the sites. The ques 
fon ^whether the Delaunay triangula- 
tion in the plane is Hamiltoman arose in 
pattern recognition problems where 
a reasonable simple curve through 
the given sites is desired. Although 
counter-examples were grjjby 
Kantabutra [1983] and Dillencourt 
[1987a, 1987b], Delaunay triangulations 
have been used with success for this 
problem. This gives evidence that this 
kraph is Hamiltonian with high proba- 
bility. Dillencourt [1987c] supports this 
thesis by proving a result concerning the 
connectivity of Delaunay tnangulation^ 
A (euclidean) traveling salesman tour 
is a minimum-length cycle spanned by n 



Figure 28. Delaunay triangulation and salesman 

tour [Dillencourt 1987]. 



sites Constructing traveling salesman 

ours is a prominent problem 
torial optimization. Pa pa dimitnou [1977 
showed its NP-completeness. The prob 
Lm becomes tractable for tours being 
cTse to optimal. It follows from the above 
'discussion that the P^^^. 

tion will ^^^^So^- 

r 2 rajs i:t^ fT ^ *. 

tour are drawn solid and dashed re- 
spectively. Rosenkrantz et al. [1977] 
howeverf observed that traversing the 
minimum spanning tree twice will pro 
duce a tour that is within a tactor oi 
2 of the optimum: Since removing , an 
edge from the traveling salesman tour 

S a mfnimum-length matching of their 
endpoS Christofides [1976] improved 
this to a factor of 1.5, with the expense of 
In CHn» ^ lo^ ») construction time due 



10 This is even true for special traveling salesman 
able as the intersects graph of disks aro 
^^SS^SS^ et al. 
1987]. 
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Figure 29. Relative neighborhood graph and 
Gabriel graph. 

to the best known matching algorithm by 
Vaidya [1988b]. 11 

2.4.3 Relative Neighborhood and Gabriel Graph 
It is worth mentioning two graphs that 
ll\n between" the mining = ng 



fortunately, no algorithm with signifi- t 
cantly subquadratic running time is ; 

T"shnilarly defined construct is the 
Gabriel graph. It contains an edge be- • 
tween two sites if the disk having that 
edee as its diameter is empty of sites. 
This concept proved useful in processing 
geographical data [Gabnel and Sokal 
1969; Matula and Sokal 1980] Howe 
U978] observed that the Gabriel graph 
consists of just those Delaunay edges that 
intersect their dual Voronoi edges Based 
on the same observation, Urquhart [1980] 
gave an 0{n log n) time algorithm for its 
construction. It is easily seen that the 
Gabriel graph is a supergraph of the rel- 
a£ve neighborhood graph. The edges of 
?he former are shown solid or dashed m 
Figure 29. 



S 1 



Tee andTe Dd^ay triangulation in 
he sense of making explicit more pro^ 
unity information than does the former 
Sless than does the latter.. Among 
i th Relative neighborhood graph 
^connects two sites provided no other 
Se is closer to both of them than thei 
interpoint distance. In Figure 29, the 
edges of the relative ^neighborhood graph 
are shown solid. Toussamt [19801 re 
ported the usefulness of this structure in 
pattern recognition, and Supowt [1983 
succeeded in constructing it in 0(n log n) 
Sme given the Delaunay trian^la -n 
Yao (personal communication, 1988) an 
nounced that 0(n) time suffices for de- 
riving the relative neighborhood graph 
Som the latter. Interestingly, this graph 
retains its linear size in higher dimen- 
sions (provided a general placement of 
he sites). This stresses the 
its rapid construction, especially m three 

j- ««i«nq since it would imply fast 
dimensions, sini-e ^ •j.\, mc Tin- 

minimum spanning tree algorithms. Un 



2.5 Clustering Point Sites 

Clustering a set of n point sites in the 
pUne melns determining partitions of 
the set that optimize some predefined 
clustering measure. This measure usu- 
ally is a function of the mterpomt 
distances of the set to be clustered. 



11 



the Delaunay triangulaticms of the sites. 



2.5. 1 Hierarchical Methods 
A single linkage clustering JHartigan 
1975 hiearchically clusters the sites as 
follows. Initially, the sites are ^nsidered 
to be clusters themselves. As long as 
there is more than one cluster the two 
closest clusters are merged. Th dl ™ 
of two clusters is defined as the mini 
mum distance between any two sites, one 
from each cluster. This process can be 
carried out in a total of 0(n log jOtuje - 
by maintaining, at each stage, the fol 
lowing Voronoi diagram: Each point ot 
the plane belongs to the region , oi .the 
closest cluster, and closeness is w th re 
sp ect to the closest site within _a cluster 
The regions of the two closest clusters 
3 always have a common edge m the 
diagram. That edge corresponds to an 
edge of the minimum spanning tree ot 
2f entire set of sites [Shamos and Hoey 
1975]. 
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Figure 30. Complete linkage Voronoi diagram. 



The model of complete linkage cluster- 
ing is similar, except that the distance of 
Z clusters is defined as the — 
distance between any pair of sites one 
from each cluster. The corresponding 
f Cnoi diagram, where closeness o a 
point in the plane is with respect to the 
furthest site in a cluster, deserves specia 
SSentSn. Its size is fi(» 2 ) for general 
clusters and 0(a) provided the convex 
hulls of the clusters are pairwise disjoint 
ffidelsbrunner et al. 1989]. Figure 30 
Sow the diagram induced by four two- 
1 clusters, indicated by bold segme n s 
For each cluster, the portion of the ! bisec 
tor of ^ two sites that lies within its 
region is drawn dashed. The diagram can 
be constructed in 0{n log n) random- 
ized time for convex-hull disjoint and 
constant-sized clusters and- m this 
setting- applies to the [ retrieval of 
expected-nearest sites [Aurenhammer 
St al 1991]. Efficient maintenance ot the 
diagram during the clustering process 
would lead to an improved clustering al- 
gorithm but remains open. 

2.5.2 Partitional Strategies 
As another example, consider the prob- 
lem of partitioning a set of n sites > into , a 
fixed number of t clusters such that the 
minimum of the (single linkage) dis- 
tances between clusters is maximized. 
The following elegant 0(n log n)-time so- 



lution was proposed by Asano et ah 
[1988] Construct the minimum spanning 
tree of the sites and remove its t -I 
longest edges. The resulting t subtrees 
LlSat Sve the desired clustering Us- 
ing similar methods, Asano et al [19»»1 

attain an equally ^^\^^ 0 
partitioning the sites into two clusters so 
that their maximum diameter is mini 
mized. This problem becomes ^consider- 

ably harder if ^"X.'J 
two clusters are sought, ^rtitions 
separable clusterings have been investi 
gSed by Dehne and Noltemeier [1985] 
and Heusinger and Noltemeier [1989] 
who exploited the fact that clusters of* 
Its being separable from .the remaining 
sites by a straight line define unbounded 
region! in the associated order-* Voronoi 
diagram. 



2 5.3 Optimum Cluster Selection 
The k-variance problem asks for select 
ins from a given set S of n sites, a 
Ster of k sites with nrijarjj 
ance the sum of squares of all mtersite 
"stances in the cluster. Aggarwal et ah 
ri989b] who cite applications in pattern 
recognition, observed the following inter- 
est fact.' If C is ^-cluster of minimal 
variance, then the region of C m the 
order-fe Voronoi diagram ot & _ is 
nonempty. Thus it suffices to examine 
The variances of sets of k sites associated 
wfth Voronoi regions. This can be done 
Tn time proportional to the number 
0(k(n - ft)), of regions by exploiting the 
SXt siis of neighbored regions differ 
in exactly two sites. Thus, the most 
time consuming step is the precomputa- 
tTn of a higher order Voronoi diagram. 
A similar approach leads tc . improve 
ments if the cluster measure to be mini 
mized is not variance but diameter. If C 
5 a ft-cluster with minimum diameter 
then C is contained in some set of 3ft -6 
sites whose region in the order-(3& - 6) 
Voronoi diagram of S is nonempty; 

^I^rTet al. [19881 claimed its NP-completeness^ 
ACM Computing Surveys, Vol. 23. No. 3, September 199! 
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see 



Aggarwal et al. [1989b] where 
7n 0(k 2 nJk log rc)-time selection 
algorithm is obtained in this way. 

Several other clustering algorithms 
that rely on Voronoi diagrams have been 
described in the literature; some of them 
are mentioned in the Introduction. 

3. SELECTED TOPICS 

A thorough understanding of the geomet- 
ric combinatorial, and topological prop- 
erties of Voronoi diagrams is crucial for 
the design of efficient construction algo- 
rithms. This third part of the survey 
presents these properties in a unified 
manner and discusses their algorithmic 
implications. We start with a description 
of the various relationships of Voronoi 
diagrams to objects in one more dimen- 
sion We continue with some topological 
properties that are particularly useful for 
the divide-and-conquer construction of a 
large class of Voronoi diagrams in the 
plane. Finally, we investigate a continu- 
ous deformation of the planar Voronoi 
diagram being well suited to construc- 
tion by the plane-sweep technique 
and extend it to generalized Voronoi 
diagrams. 

3 1 Geometry of Voronoi Diagrams: Their 
Relationship to Higher Dimensional 
Objects 

Voronoi diagrams are intimately related 
to several central structures in discrete 
geometry. This section exhibits these re- 
lationships in a comprehensive manner 
and discusses some of their geometric 
and algorithmic consequences. Although 
our main interest is in the classical type 
of Voronoi diagram, it is advantageous to 
base the discussion on the more general 
concept of power diagram since the geo- 
metric correspondences to be described 
extend to that type in a natural way: We 
shall refer to d dimensions m order 
to point out the general validity of the 
results. 

3.1.1 Convex Polyhedral Surfaces 
Consider a set S of n point sites in R d 
such that each site p e S is assigned an 




Figure 31 . Power line of two circles. 



individual real number w(p) called its 
weight. The power function of a point 
xeR d with respect to p is specified by 

p0 w(*,p) = (x-p) T (x-p) - w{p). 

For w(p) > 0, one may think of the 
weighted site p as a sphere in K with 

radius v^(pT around P; f0 \ a ^ p rla 

outside this sphere, pow(x, p) > 0, and 
J^UTp) expresses the distance of x 
to the touching point of a line tangent to 
the sphere and through x. In the plane, 
the locus of equal power with respect to 
two weighted sites p and q is commonly 
called the power line of the correspond- 
ing circles. If the circles are nondisjoint, 
their power line passes through their 
points of intersection; see Figure 31. In 
general dimensions, the locus of equal 
power is a hyperplane m 2? , called the 
chordale of p and q or chor(p, q) , tor 
short. chor(p, q) is orthogonal to the 
straight line connecting p and q but 
does not necessarily separate p from q. 
For w( P ) = w{q), chor(p, q) degenerates 
to the symmetry hyperplane of P and q. 
Let h(p,q) denote the closed halfspace 
bounded by chor( p, q) and containing the 
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Figure 32. Spatial interpretation of chordale [Aurenhammer 1987a]. 



ts of less power with respect 
power cell of p is given by 

cell(p) = f] h(p, q). 

qeS-{p} 

lalogy to classical Voronoi regions, 
lower diagram of S, PD(S) is the 
;x polyhedral complex defined in R d 
ese cells; see Figure 4 for the case 
. PD(S) coincides with the classical 
loi diagram of S in the case of 
ly weighted sites. Note that power 
may be empty if general weights 
;ed. 

technical reasons let us identify 
ith the linear subspace of R d+1 
;onal to the (d 4- l) st coordinate 
: d+v The key observation for the 
ice of objects in R d+l that are re- 
o power diagrams is that the power 
>n pow(x, p) can be expressed by 



to p. the hyperplane 

k{p)- x d + i = 2%t p ~ p t p + w (p) 

in R d+1 [Aurenhammer 1987a]. This fact 
is best explained by reference to Figure 
32. Let p and q be two sites in it! 2 
whose weights are indicated by circles. 
Projecting the circles onto the paraboloid 
of revolution x 3 = xf + x\ defines two 
planes ir(p) and r(q) that intersect the 
paraboloid in these projections. It is an 
easy analytical exercise to prove that the 
line 7r(p) Pi ir(q) projects orthogonally to 
the chordale chor(p, q) in the x 1 x 2 -plane. 
This correspondence holds for arbitrary 
dimensions. As an immediate conse- 
quence, PD(S) is the orthogonal projec- 
tion of the boundary of the polyhedron 
that comes from intersecting the half- 
spaces of R d+1 above x(p), for all peS. 
In fact, x is a bijective mapping between 




^C;;;'^^^'^'^\'.'■ v ^,;^l,'■^,V■.^^ ; . 1 >\ , ^^'•.■ , . .''\\ '" 
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weighted sites in R d and hyperplanes in 
Thus, for any intersection of upper 
half-spaces there exists a corresponding 
poweT diagram in one dimension lower 

a 1htVene S ral result has far-reaching 
impYicatfons. Most important, it shows 
that power diagrams are m a geometric 
and combinatorial sense equivalent to 
unbounded convex polyhedra or more 
Drecisely, to their boundaries, that is, to 

hedra are well-understood objects 
in discrete geometry, hence so are power 
diagrams. Exact upper and lower bounds 
on L numbers of their faces of various 
dLensions are known. ^ particular 
PD(S) in R d realizes at most n ceus, jj 
./-dimensional faces {l*j*d- D, and 
f 0 - 1 vertices, for 

n - d + i - 2 

i 











d 


a = 


2 




2 . 



This follows from the so « - 
bound theorem [Brondsted 1983]. The 
numbers fj are 0(n ^'^ ) for 0 < j* 

d As 1 for the facets of a convex polyhe- 
dron the power cells are convex but pos- 
S' unbounded polyhedra. Not every 
half-space within a given set needs to 
contribute to a facet, so power eel s may 
be empty or degenerate, d + 1 or more 
hyperplanes in intersect in a com- 

mon point (unless parallelism occurs ) so 
«t lpast d + 1 power cells meet at each 
vertex of PD(S). In particular, there are 
no vertices if the cardinality n of S does 
not exceed d. These and many other 
properties of PD{S) can be read off con- 
veniently from its (d + l)-dimensional 
embedding. 



3.1.2 Hyperplane Arrangements 
In order to exploit the situation fully, 
we turn our attention to diagrams of 
higher order [Edelsbrunner 1987), 
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Aurenhammer [1987a]. To this end, the 
concept of a power cell is generalized to 
more than one site. Let T be a subset of 
k sites in S. The power cell of T i s 
defined as 

cell(T)= H HP-<7)- 

peT,qeS-T 

The complex of all (nonempty) power cells 
arising from the (j) subsets of S with 
fixed cardinality k is known as the 
order-k power diagram of b k-FD(i>). 
Clearly l-PD(S) = PD{S) holds Note 
that k-PD(S) is just the order-fe Voronoi 
diagram of S provided all sites in S are 
equally weighted, (n- iyPD{S) .is also 
called the furies* site power diagram 

f s 

Now let Z denote the intersection of 
the half-spaces below the hyperplanes 
x( p) for all p 6 T and above the hyper- ; 
planes HQ) for all _,eS - T .Due to 

ffiSSSon^a^. connection 
between order-fe power diagrams and so- 
called hyperplane arrangements, lne flat- 
ter are cell complexes that arise from 
dissecting the space by a finite number of 
hypSlanes. Each cell in an arrange- 
ment is the intersection of some half- 
snaces below or above its defining n hy- 
S pSSnes. If we let a consist of aU 

Sells lying below k and a bove r, i - k hy- 
npmlanes we can conclude that the « 
Lvel Scedby {*(p)| peS> projects to 
the cells of k - PD(S) for each fc be- 
tween 1 and n - 1. This shows that each 
Sraiement in corresponds to a 

SSTfamilj of higher order power 
diagrams in R d and vice versa. The cor 
Smdence between order-1 power d . 
Krams and convex polyhedral surfaces 
discussed earlier is a special case of this 

" F$ure 33 depicts a cell (left) being the 
infeSection of five upper and three owe 
half-spaces, that is, a three-level cell. Its 
orthogonal projection onto the 
is an order-3 power cell among eight sites 
right) whose boundary is indicated by 
bold ines. Note that the power cell is 
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Figure 33. Projecting an arrangement cell [Aurenhammer 1987a]. 





carved into subcells by the furthest site 
power diagram of the three emphasized 
sites and by the power diagram of the 
remaining five sites. 

The correspondence to arrangements 
straightforwardly yields an upper bound 
on the overall size of all higher order 
power diagrams of a fixed set of n sites: 
the maximal number of faces of ^ar- 
rangement of n hyperplanes in R 
This number is precisely 



d-j + 1 



or faces of dimension j, which is in 
l(n d + l ) for 0 <j < d + 1 [Alexanderson 
ind Wetzel 1978]. Edelsbrunner et al. 
1986a] developed an algorithm for con- 
tructing such arrangements in 0(n ) 
ime. Consequently, the complete family 
f higher order power (or Voronoi) dia- 
rams of a given set of sites can be 
Dmputed in optimal time and space. 
Idelsbrunner [1986] showed how to con- 
duct the fe-level of an arrangement 
i R 3 , and thus a single order-fc 
Dwer (or Voronoi) diagram in the plane, 
: a cost of 0(y/n log n) per edge, 
ecently, Mulmuley [1991] proposed a 
tndomized algorithm for constructing 
vels of order 1 to k in d > 4 dimen- 
ons in worst-case optimal time. 



We state another interesting conse- 
quence. Let x'(p) denote the image of 
reflection of the hyperplane ir(p) through 
R d . The cells in the fc-level of {ir'(P) I P e 
S\ correspond bijectively to the cells m 
the (n - *)-level of {r{p)\ peS} by re- 
versing "upper" and "lower" for their 
supporting half-spaces. Clearly, any two 
corresponding cells yield the same or- 
thogonal projection onto R . That is to 
say, they yield the same power cell that 
is simultaneously of order-fe and order- 
(n- k) for fixed k. This shows that any 
order-fe power diagram for n sites can be 
intepreted as some order-(n - fe) power 
diagram. In particular, the furthest site 
power diagram, of S is the order-1 power 
diagram of some set of sites (distinct from 
S, in general). We point out that this is 
not true for the subclass of classical 
Voronoi diagrams, as can be seen from 
the maximal size of their closest site and 
furthest site counterparts that differ even 
in R 2 \ see the Introduction and Section 
2.2.1. 

Finally, let us mention a result on 
power cells in R 2 . Let T C S and \ T\ = 
k By definition, cell(T) is the intersec- 
tion of k(n - k) half planes and thus 
might have that many edges. On the 
other hand, cell(T) is the projection of a 
polyhedron in R 3 that has at most n 
facets and thus only O(n) edges by 
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in 0{n log n) time by intersecting 
half-spaces of R . 

3 7.3 Convex Hu//s 

The conuex fcuH of a finite point set is 
lhe convex / half-spaces con- 

^ ^Xset N^t we stud? the dual- 
Sefl) Power diagrams ? nd 
clvex hulls ^Aurenhammer 198 a]. 
Duality is usual y defined to 

vex ^ ei \l h ll\J R ^ Z and Z' 
r rdlo^^anf thefe is a bijective 

faces of Z and lo J 
faces of Z' such that / £ g, } - t , ta \ a 
Sees f and g of Z if and only if Hf IS 
fr A Since power diagrams in B are 
Vl/ )• ainct; rav r.oivhedral surfaces 
projections of convex P 0 ^ 6 "?' ies 
in R d+ \ the notion of duality carries 
over in a natural manner. 
° We now consider, for each site peS 
with weight u,(p), a point m K ■ 




Figure 



34. Polarity between points and lines. 



A(P) = ( 



p T p - w{p) 



From the formula expressing the ^hyper^ 

hi the geometric interpretation dis- 

h f J i •„ Fimire 34 The hyperplanes 

F^htrtheTraboloid at its points of 
touching tne pai <* A( ^ 

intersection with ir(p) au cou^ 

Recall that the inters ec tion, Z o f the 

h alf-spac-es above .(p) for all PjS, P™ 

-A kdbt " P "e- 
planes defined by a. o"«-f > f , i „ n 

=lrv PS the relative position between 

Sts and hyperplanes, there is a hyper- 
points ana nyp F lt and 

nlane with A(Pi). •■•> ayPm '. :„„ 
t{q) above it, for all remaining sites 



ae S This, however, means that 
\ A( P ) span a (d - j)- 

Sat aS visible tarn £ 

laceb ui 1.1 _ O0 j n otner 

point on the 3C d+ i' axls ,. f 7 

words we have shown the duality oi Z, 
IS thus of PD(S), to that lower convex 

hU By P the same reasoning, the duality of 
the furthest site power diagram of S and 
lug f ,i convex null oi 

% an be established. This 

^I'a^^** link between 
SSst site and furthest 
grams and, particularly, classical 
SagramB. Most importan hese d^ 

d u- L hoth of which are worst-case 
achieved, ootn oi wu. „ ri9771 
optimal; see Preparata and Hong ; LIS ' ' ^ 

i y P /nTunpfr part of the convex hull 
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triangulation whose edges are orthogo- 
nal to the facets of a given cell complex 
in R d i s a necessary and sufficient condi- 
tion for that cell complex to be a power 
diagram or, equivalently, to be the pro- 
jection of a convex polyhedral surface 
[Aurenhammer 1987b, 1987c]. If the 
given cell complex is simple (i.e., exactly 
j+l cells meet at each vertex), a projec- 
tion surface can be reconstructed, if it 
exists, in time proportional to the num- 
ber of facets of that cell complex. It is 
veil known that such a surface always 
ixists except for d = 2. This implies that 
xty simple complex of n cells in d >: 3 
imensions, which might have a large 
umber of lower dimensional faces, can 
e specified by storing only n half-spaces 
7 R d+1 

Surprisingly, even higher order power 
agrams are dual to some convex hull in 
d+1 [Aurenhammer 1990a]. For any 
bset T of S with k sites let A(T) 
note the point in R d + 1 whose ortho- 
nal projection onto R d and whose 
+ l) st coordinate are 

£p and 5^ [p T P ~ w (p)}> 

peT peT 

pectively. This definition conforms to 
case k - 1 since A(T) = A(p) when 
= { p} . Using arguments similar as for 
1, one may verify that the lower part 
he convex hull of the point-set 

S k = {A{T)\TCS,\T\ = k] 

lual to k-PD(S). As one corollary, 
very lower part has to be dual to the 
rsection of some corresponding upper 
spaces and hence to an order- 1 power 
ram. To be more precise, k-PD(S) 
:ides with PD(S' k \ for S' k being the 
•rdingly weighted) orthogonal projec- 
of S k onto R d . This implies that the 
of power diagrams is closed under 
'-k modification— another result that 
t true for the more restrictive class 
issical Voronoi diagrams. From the 
ithmic viewpoint, k-PD(S) becomes 
ructable via convex hull algo- 
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rithms. For example, the planar Voronoi 
diagram of order-& can be obtained from 
its order-(& — 1) predecessor in 
0(kn\og n) time by simply computing a 
convex hull in R 3 . Insertions and dele- 
tions of sites in an order-& Voronoi dia- 
gram also amount to computing certain 
convex hulls; these operations even can 
be handled on-line in an efficient manner 
[Aurenhammer and Schwarzkopf 1991]. 

3.1.4 k-Sets 

Let Q be a set of n points in R d+1 . A 
k-set of Q is a subset of k points of Q 
that can be separated from the remain- 
ing n — k points by some hyperplane of 
R , The discussion in the foregoing 
section reveals a correspondence between 
fc-sets and order-& power diagrams. 

By polarity, points in Q can be trans- 
formed into hyperplanes; each point q e Q 
bijectively corresponds to a weighted site 
p in R d such that A(p) - q and thus to 
the hyperplane 7r(p) in R +1 . Now, let 
M be a k-set of Q, and consider a hyper- 
plane h that separates M from Q — M. 
Without loss of generality, M lies below 
h. Equivalently, the pole of h lies below 
all polar hyperplanes of the points in M 
and above all polar hyperplanes of the 
points in Q — M. Hence the intersection 
of the corresponding lower and upper 
half-spaces of R d+1 bounded by these 
polar hyperplanes is nonempty. This in- 
tersection, being a cell in the &-level of a 
hyperplane arrangement, projects to a 
cell of the order-& power diagram of a set 
of weighted sites, { p e R d \ A(p) e Q} . 

We have obtained a one-to-one corre- 
spondence between k-sets in R d+1 and 
power cells in R d . This allows us to ap- 
ply known results on the number of k-sets 
to the analysis of the size of order-/e power 
diagrams. Let f d {k, n) denote the maxi- 
mum number of &-sets of any set of n 
points in R d . Known asymptotic values 
of this function are f 2 (k, n) = Q(n log k) 
and f 2 (k t n) = 0(ny/k) [Edelsbrunner 
and Welzl 1985], f 3 (k, n) = Q(nk log k) 
and f 3 (k, n) = 0(n 8/B log 5/3 n) 
[Edelsbrunner et al. 1986; Aronov et al. 
1990, respectively]. Not surprisingly, a 
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better understood quantity is 



g d {k,n)= ZW>") 



Chazelle and Preparata [1986] showed 
Chazeue an ciarkson [1987] 

¥ \5 ~a (k n) = 0(n^k>), for any 
showed g 3 l«> ^\ i „^ - 0(n d ) 

[1989] succeeded in proving g d {k, n) 
ruTl*' 2 * fe rd/21 ) as n/fe-°° f° r tixed 
^Determining exact or at least sjyjg 
totically tight bounds on f£K ^ 
remains an important open problem. 



higher dimensions. Below we will show 
that a certain rather general type of 

T7 ; Jiatrrnm in B can be embed- 
Voronoi diagram in ^ , 

in to a power diagram in K . lo this 
end we identify B * with the linear 
inace x„ = 0 of B d . Let p be some 
noS s e ?n B- 1 and consider an arbi- 
trary distance function flx, p) for p For 
some strictly increasing function F on 
i ! we define the com of p with respect 
to F as 




3.1.5 Related Diagrams 



TW central role of power diagrams within 
^ context of Voronoi diagrams becomes 
the contexx in observations 

1 Li sites For the three chordales 
T g chor(p, q) 0 chor(q, r) C 

they d f n n % e JS^ h 9 0 \ds. It is not dif- 

Sdales in B d and intersections of hy- 
chordales m xv her hand) 

perplanes m B ^ fey the 

^^r^Sor Ly VoUoi diagram with 
'ThedS regions since the regions 

!™ wioS separators are hyperplanes 
aSarn TMs resu>« -PP J- 
to diagrams defined by tne & 
quadratic-form distance 



xsR d -\x d -F{f(x,p))]. 

A Voronoi diagram for point sites in 
pi-i un der the distance function f is 
fermed transformable if there exists an 
F that forces the affine hull a(p, <Z), of 
conejp) n cone F (q) to be a hyperplane 
C0Ue Id P L » nv distinct sites p and q. 



"j i, ^ 



pone^p) n cone F (,q; "° - r . 
in B d for any distinct sites p and q. 
transformable diagrams represent pro- 
ved sections of iT^^f< 
* f and only if its P-jection ^, 
£to cone.(p) lies in a fixed (open) ha£ 

rt >,/n bounded by a(P, 9)- AS 
space, h{p,q), ~u r ion 

a consequence, x belongs* ^ 
exactly if 



x'eZ= 0 MP.9)' 

qeS-{p) 

f«r <? being the underlying set of sites in 
- She d-dimensional polygon 

S is a power cell since a(p, q) n «(«, r) 
r !( /r) holds so that these affme hulls 
elAmterpreteda,^^- 
elusion, the boundary of psregiw 

orthogonal P^^"; ce ll in 
cone F (p) intersected with a powe 

E Transformable Voronoi 
for the additively weighted distance, 




■VV.°-' 



* 



with M a nonsingular and symmetric 
7Z dVmatrix and thus to ^PjJ 
i«r distance functions considered in the 
UWure Note that Q equals the power 
literature, i identity matrix. 

fU One° is tempte? o believe that power 
diagrams Related solely to objects m 
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and the multiplicatively weighted dis- 

tance, 

«KpI wit hF(m)=2m 2 
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Figure 35. Cones for two distance functions. 



The former gives rise to hyperbolically 
bounded and star-shaped regions (the 
Johnson-Mehl model, Figure 6), whereas 
the latter induces spherically bounded 
regions (the Apollonius model, Fig- 
ure 12). Its regions are disconnected, in 
general, which already reflects the major 
complexity of the latter type. 

Figure 35 displays cones for the addi- 
tively weighted distance (above) and the 
multiplicatively weighted distance (be- 
low). They are cones of revolution in the 
former case— intersecting in a hyper- 
bola that projects down to a hyperbolic 



separator— and paraboloids of revolution 
in the latter case— intersecting in an el- 
lipsis that projects down to a circular 
separator. In both cases, the affme hull 
of cone intersections is a plane. 

Upper bounds on the size of both types 
of diagram in d dimensions can be de- 
rived from the maximum size of a power 
diagram in R d+1 . Concerning the algo- 
rithmic aspects of this relationship, the 
only known algorithms for constructing 
these diagrams in d > 3 dimensions are 
obtained. The Apollonius model in the 
plane can be computed in 0(n 2 ) time and 
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generalizes nicely to order k. 

3.1.6 Upper Envelopes 

_ " c f be piecewise linear d- 
Let ■/!>.••> h r nrf T n the follow- 
variate functus on , R J- & 
ing> we will not dytm^rn Rd+1 

SSSi^poSli nfaxLum, * of 

these functions: 

U( x ) = max 

rf is called the apper f™ eto pe of 
? f We have already observed that 

'i' " * *. r n^Tmlnnpq are related to 
certain upper J™ set S of 

power diagrams. F or ^ any 
weighted sites in R . the upp 

1™ o^eronoi ^^S^ 
p,L and consider a P«Jf Ce 

cluster C, is defined as 

^ X)Cl ) = max{5(x I p)|peC 1 }. 



h is also called the ^f^^. 
tance; compare Section SUU. Each ^clus 

^ ? ryJSS^O ^B 8 , being 
polyhedral surface cap^J that 

the boundary of th e ] o aces 
comes from intersecting 
below *(p) for alj P e f ^ "^t 

sects cap(.Oi) m d y u m>er en- 

r/-' \ TVii<? implies that tne uppei 
cap(C,-). Ihis impi pro jects to 

velope of ^^r-l^^Lcei by 
^"rnd ? ffan illustration of 
Se Voronoi diagram for two-site clusters 

See T rob r t e ain° another .-niple^sume 
that visibility of the sites is constrained 

order-fe Voronoi diagram ot £>- 



■ J,,,., ty a segment on a line 
t0 convex hull of S. A site 

D Ts and a pS « are called «i.iWe if 
fhe line segment joining « and p inter- 
sects W The peeper's Voronoi dtagram 
f I and W assigns each point in the 
\ /to the region of the closest site 
Se from it lee Figure 36 where the 
window^ shown between two bold line 
window diagram is composed of 

Spen^uS bisects of sites (dashed) 
and of rays through sites and window 

e 1SStf pmay be associated with an 
unbounded convex polygon plat e(p) in 
$ being the projection onto *(p) of the 
' TiflJ, visible from p. Again it 
set of all points J 15 ^ f u into the 
is easy to see that a pointy ^ 

+w hps above plate(q) to r au * c u 
We conclude that the upper enve- 

0 P p } e of Plate(p) • .^VoTs 

Upper env v tudied ge0 met- 

functions on K are w 
ric objects. Pa* and Sh a nr U ^ J 

"t^is^loS number of triangles 
where n is tne ^ linear pieces 

needed in Phoning the hnea p 

of these ^p^ma^s function, 
the inverse of Ac*wm . t that run 
Construction ^ gonthms exu * 
in time proportional to their 
size; see Edelsbrunne et _aL U^l 
also give a lower bound ol ° _ 

„f « cluster Voronoi diagram. Ke 
size of a cluster enve l 0 pes by 

cent results on ^ u T that the 
Huttenlocher e ; J^f ^ t H Lai of » 

di f ^nTs fsi e Tol^n)). The upper 
»tes has a size ot V ^ pe 

envelope ge size 0 f 

Voronoi diagram has a wo 
CXu 2 ) [Aurenhammer and StocKi i» 
xt * iCni the superlinear size ot a 
Note that .^ h D e lie S s P the existence of 
diagram implies 

disconnected regions. 





3.2 Topology of Planar and lt8 

Divide-and-Conquer Construction d 

Variants 
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Figure 36. Peeper's Voronoi diagram. 



ght into connection with geometric 
ts having nice algorithmic proper - 
like convex hulls or hyperplane ar- 
^ments. So there is need for finding 
res common to such "ungeometric" 
ams. We do this by extracting topo- 
il properties, being guided by their 
mce to the divide-and-conquer con- 
ion of Voronoi diagrams. Attention 
tricted to the plane as the instance 
3t applications and easiest analysis. 
! concept of diagram we are inter- 
in here may be characterized as 
s. The n sites are points in R 2 . 
istance function varies by different 
5 of its "circles" that may depend 
al properties of R 2 . It thus may 
Dr the individual sites. 

Convex Distance Functions 

portant type fitting into the con- 
bove is generated by the convex 



distance function or the Minkowski dis- 
tance [Chew and Drysdale 1985], Let C 
denote the boundary of some compact and 
convex subset of R 2 with the origin o in 
its interior. The distance with respect to 
C of o to some point x e R 2 is given by 



d c (x, o) = 



d(x ) o) 
d(x\ o) 



for x r being the point of intersection of C 
and the ray from o to x. Clearly d c can 
be defined with respect to any site q of a 
given set S by translating C so that q 
occupies the former position of o. Ob- 
serve that d c (x, q) depends on the eu- 
clidean distance 5(x, q) as well as on the 
direction of the ray from q to x. Convex 
distance functions include, among others, 
the general L p -metric 



& P {x,q) = v / l<h-*ir+ 1<7 2 - 
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and x = 




whose most important instances are for 
„ - 1 (the Manhattan metric; C is a 
^arerotated by 46'). : P = 2 e«- 
Sean dirfonce; C -cle) and P 

^anKS SS^uS £ is a met- 
parallel square^ illt . sym metric 

nc only ii ° l * , , mieht 
with respect to q since dc(9,r ^ign 
differ from d c (r, 9 ) otherwise. However, 
the triangle inequality 

d c {q,r) + d c (r,s)>d c {q,s) 

can easily be shown to hold even if C is 

? and region reg(q) of q; compare the 
iuctfon. Clearly sep(, r need n 
Hp a straight line as is the case 101 
idtdefn fiance. A». even more «n ; 

ring to the lexicographical order (sym 
hSized bv <) of the sites m S [Klein 
and Wood 1988]. The dominance of q 
over r is redefined as 





Figure 37. Redefining a separator. 



Dom(q, r) 



dom(q.r) -sep(g.r), 
dom(q, r), 



if q < r 
otherwise 



The new separator, Sep(q, r) , is the in- 
tellection of the boundaries of Dorn(q, r) 
STSSn(r, 9 ), which is obviously one- 
dimensional. Figure 37 shows a two 
dimensional separator m the L.-Jj™ 
(left) and its redefined version (right)- 
For a Voronoi diagram in the L^metric, 
see figure %. The region of a site , now 
can be rewritten as 

Reg(q)= H Dom(q.r). 

reS- {q} 



That is, each point in R 2 belongs to the 
region of the lexicographically lea 
Z am on2 its closest sites in S. This 
ensured that the interiors of the regions 
^over the plane up to a set of Lebesgue 

m De S fined e L° this way, the Voronoi dia- 
trn V ( ^ for a convex distance func- 
ti'on /fas he nice property that its 
reins' are star shaped; Reg(q) contains 
S^ght-line segment between 9 and 
* for eadi *eReg(<7). To prove this by 
contradiction, assume the existence rf, 
point , on this -gment w lt h^c(^) > 

/+ d Jy, r) By the triangle inequality, 
ll^tir sum iJ greater than or equal 
d r (x, r) so that d c {x, q) > d c {x, r), 
;= Y rfRe£(o), follows. 

plies their simple connectedness . Clear Y 

regions Since ^^^J^J^SS 
ST^^iS « the size of 
V C (S). 



I 



' ' . v j 





3 2 2 Divide-and-Conquer Construction 

can bT constructed efficiently by a 

their solutions recursively (unless they 
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_„ v ery small and can be solved by 
Jrivial methods), and finally combines the 
partial solutions to the g obal one The 
V type of diagram is well suited to at- 
tack by this strategy. The underlying set 
« I of sites is partitioned into two subsets 
c and S 2 of nearly equal cardinal- 
ity VASJ and V C (S 2 ) are computed 
recursively and are then merged to 

V ' The merge step deserves the main at- 
tention since it actually constructs the 
diagram. We will discuss this process in 
some detail and under the simplifying 
but computationally unrestrictive as- 
sumption that S 1 and S 2 are separable 
L a vertical line. Since V C (S 1 ) as well as 
V„(S,) defines a partition of R , any 
point x e R 2 falls into the region Re gl (g) 
nf VAS,) for some site qeS 1 and into 
the region Reg 2 (r) of V C (S 2 ) for some 
site reS 2 . Now observe that xeReg(q) 
of VAS) if d c (x, q) < d c (x, r), and x e 
Reg(r) of VAS) otherwise. So we have to 
cut off some part from Re gl (q) and from 
Reg,(r) by means of Sep(q, r) m order to 
obtain Reg( 9 ) and Reg(r). Carrying out 
this task for all relevant pairs of sites 
constitutes the merge process. The union 
of the newly integrated pieces of separa- 
tors will be called the merge chain, 

M(S V S 2 )- . , ,. r 

The efficiency of an implementation ot 

the merge process critically depends on 
the topological properties of the merge 
chain. We first show that M(S V S 2 ) be- 
haves well if the euclidean distance 8 is 
taken for d c [Shamos and Hoey 1975]. 
As a matter of fact, each horizontal line 
intersects Af(S lP S 2 ) in exactly one point 
in this case, that is to say, M(S V i> 2 ) is 
(vertically) monotone. M(S V S 2 ) is a 
polygonal line composed of edges of the 
classical diagram V(S) that are portions 
of perpendicular bisectors of sites in S x 
and S 2 , respectively. Figure 17 gives an 
illustration. If there were a horizontal 
line intersecting the merge chain in two 
or more points, some site q e S 1 would be 
forced to have a larger x r coordinate than 
a site reS 2 . This contradicts the as- 
sumption that S x and S 2 are separable 
by a vertical line; compare Figure 38. 




Figure 38. Impossible merge chain. 



The monotonicity of M(S U S a ) implies 
its connectedness and unboundedness. 
Hence the merge chain may be con- 
structed edge by edge, starting with an 
unbounded one. As is mentioned m the 
Introduction, unbounded edges separate 
sites on the convex hull of S. So it suf- 
fices to determine a site m S x and a site 
in So that are neighbors on the boundary 
of the convex hull of S. This can be 
accomplished in O(n) time by finding a 
line tangent to the convex hulls of re- 
spective S x and S 2 . Each further edge of 
M(S lt S 2 ) is found by tracing the bound- 
aries of the current regions of V&O ana 
V(S,) in an appropriate direction until 
they intersect the actual separator. Easy 
counting arguments show that only U{n) 
edges are processed in total: Each edge 
traced is either discarded or shortened to 
a new edge of V(S). This implies an 
overall runtime of O(n) for the merge 
process. By the recurrence relation 
Tin) = 2T(n/2) + O(n) that results 
from using divide and conquer, the 
time complexity T(n) for computing a 
classical Voronoi diagram is 0{n log n). 

Let us now come back to the distance 
d r for general convex shapes C [Chew 
and Drysdale 1985]. As in the euclidean 
case, M(S V S 2 ) is composed of those parts 
of the boundaries of regions ot v c (t>) 
that separate sites in S 1 from sites in b 2 . 
Thus M{S V S 2 ) is one dimensional. Po- 
tentially, this merge chain may consist ot 
cyclic and acyclic (but then necessarily 
unbounded) topological curves m K . 
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Examples show that M{S V S 2 ) is in fact 
Connected in general. Its connected 
^ponents ^^^gj^ 
TZt^t aslSV existence 
of some cycle M and deducing a contra- 
d ctX Without loss of generahty M 
encircles a subset T of sites n &!■ 
Sv then M does not encircle any 
site in S 2 . Furthermore, there exists a 
Eftmost point xeM with seSep(r,s) 
for some r,seS 2 . Since the regions of 
VMS) are star shaped, the three open 
line segments from « to J. 
o e T, respectively, do not intersect m . 
Thus as r and s lie to the right of q, x 
has to lie to the left of q. This, however, 
SLadicts the convey of the ^ 

go5K: , aXlfnce the last three 
S are equidistant from * with re- 
P +. +~ W Figure 39 illustrates this 
Suation;^ 2?C « shown bold and 

d t££E££«* the acyclic merge 
chain M(S 1; S 2 ), one needs to detect some 
Sounded 'edge of each of ts connected 
components. The previously described 
procure of tracing Wto-*- 

^rL^^ »«** that 

works successfully in the euclidean case 
Tails to be correct for d c . To ge an alter- 
Stive method, we observe that each un- 
bounded edge has to be contained m the 
Ssection g of two unbounded regions of 
VJSO and V C (S 2 ), respectively. By a 
circular scan through these regions, each 
relevant pair Re gl («) and : Reg , r) can be 
found and intersected with Sep( 9 , r) in a 
S3 of O(n) time. If Keg^ . H i Reg^ 
n Sep(Q, r) extends to infinity, then this 
edge Tsteken as the starting edge for a 
connected component of M(by 

The linear-time behavior of the merge 
process thus can be maintained for gen- 
eral convex distance functions d c - Con- 
sequently, V C (S) requires 0(n \og ri) 
time and 0(») space for construction^ Of 
course the complexity analysis is based 
on Se assumption that d c is computa- 
tionally simple. That is, d c (q «) canbe 
calculated in 0(1) time for all q e S and 
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Figure 39. Cycles do not occur. 



x e R 2 , and, in addition, ^ separators 
can be computed and intersected m 0 1 
time. This may not be the case for 
complex shapes C, however. 

3.2.3 Nice Metrics 

We have seen that divide and conquer is 
^powerful approach to the construe ;ion 
of generalized Voronoi diagrams It also 
annlies well to sites more general than 
noFnts see Section 1.3.3. We do not pur- 
^e these modifications here however 
See they do not fit into the present 
concept of a diagram. Instead, we at 
tempt to characterize metrics in R that 
preserve the two main Property of a 
Voronoi diagram exploited in the pre 
ceeding discussion [Klein 1989; Klein and 
Wood 1988]: the connectedness of regions 
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that implies an 0(n) size and the ab- 
sence of cycles in the merge curve that 
allows an 0(nlog n) time construction. 

Let d be any metric in R 2 that in- 
duces the euclidean topology. A curve 7 
in R 2 is termed d-straight if, for any 
three consecutive points x, y, z on 7, the 
equality 

d(x, y) + d(y, z) = d(x, 2) 

iolds. Clearly 5-straight curves are 
itraight-line segments, for 5 being the 
juclidean metric. By arguments similar 
0 that used for convex distance functions 
namely using the triangle inequality), it 
an be verified that any Voronoi region 
'eg(p) f° r d is d-star-shaped: 7 C Reg(p) 
olds for the d-straight curve 7 from p 
>any xeUeg(p). The d-star-shapedness 
? /Reg(p) does not necessarily imply its 
jnnectedness since d-straight curves do 
)t always exist. If, however, for any two 
>ints x, z e R 2 there is some y distinct 
5m x, z with d(x, y) + d(y, z) = 
x, z) y then the existence of 7 is guar- 
teed. Hence the Voronoi regions under 
□se metrics are connected. 
Posing an additional restriction on the 
;tric d makes the Voronoi diagram 
m more well behaved: Assume that, 
any point meR 2 and for any positive 
mber r, the generalized disk 

D = {xeR 2 \ d(x t m) < r) 

h center m and radius r under the 
;ric d is a simply connected set. For 
d r, the shape of D may vary with 
position of m provided D remains 
ply connected. Symmetric convex dis- 
:e functions d c (i.e., with point- 
metric shape C) constitute a special 
of translation-invariant disks. And, 
ed, the proof that merge chains are 
lie can be extended from convex dis- 
e functions to d. 

underline the usefulness of this con- 
let us give some examples for d. Let 
; any point in i? 2 . The Moscow met- 
uth respect to K induces d-straight 
is composed of pieces straight to- 
K or radially around K of minimal 
euclidean length. Figure 40 depicts 
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a d-straight curve (dashed) between two 
points x and y and a Voronoi diagram 
(solid) for eight sites in this metric. For 
the geodesic distance among polygonal 
obstacles, d-straight curves are the mini- 
mum euclidean length connections avoid- 
ing any edge of the polygons (compare 
Figure 26). Also of interest are composite 
metrics where, for instance, distances are 
measured in the L r metric in some por- 
tion of R 2 and in the L 2 - (i.e., euclidean) 
metric in the complement. Here, d- 
straight curves are composed of those in 
L x and L 2 such that curve length is min- 
imized. In all three examples, the respec- 
tive disks are simply connected. As a 
consequence, the resulting Voronoi dia- 
grams—except for the geodesic distance 
function that fails to be computationally 
simple— can be computed in optimal 
0(n log n) time and O(n) space by divide 
and conquer. Thereby, the underlying set 
of point sites has to be divided into sub- 
sets by curves whose intersections with 
any disk for d is either connected or 
empty. Consult Klein [1989] for further 
details. We pose it as an open question to 
find equivalent characterizations for 
metrics whose Voronoi diagrams are 
constructable by divide and conquer in 
optimal time. 

3.3 Deformation of the Voronoi Diagram: 
Plane-Sweep Technique 

The variety of techniques that can be 
used for constructing Voronoi diagrams 
is possibly as fascinating as the diagrams 
themselves. The construction of the 
classical planar Voronoi diagram by 
means of divide and conquer or via 
three-dimensional convex hulls are 
well-investigated algorithmic problems. 
Although the achieved complexity bounds 
are worst-case optimal, researchers con- 
tinued considering alternative methods 
of construction. Aside from the incre- 
mental insertion strategy described in 
Sections 1.3.1 and 1.3.2, the plane-sweep 
technique to be discussed now stands out 
by its conceptual and computational sim- 
plicity. It extends nicely to various 
generalized types of Voronoi diagrams. 
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Figure 40. Voronoi diagram in 



the Moscow metric [Klein 19891. 



3 3 1 Deformation 

le^oU &st we introduce a contmu 
lus deformation of V(S) [Fortune 198U 
Sen we describe the plane-sweep , tecb- 
Sqle and its use for computing this 

de T°Xia^ the mechanism of deforma- 
tion the euclidean distance function with 
respect to the sites is interpreted ^ the 
Swing way. With each site peS, * % = 
cone 
cone(p) 

= {( x ,z)eR 3 \xeR 2 ,z = ?>{x>P)'l 

is associa^d. Con f ( P ) is 
rected cone m K wiui 
revolution, with apex j. ^ ™ be 
interior angle ot R 2 We 

viewed as bivariate ^°™°** e( Z I 
define the Zou^er e^eZope of {cone(p)l 



n S\ as the pointwise minimum of these 
p e 5} as uie f h surface 

functions or, equivalency , 

«^ rt f that Dortion ot eacn cone 

ssr ^ 

nroiect the envelope onto R m a miter 

in positive jc 2 -airecuuii. v */o\ n r 

e +v,p desired deformation V (&) ot 
g vS To be more precise, the deforma- 
2n* maps each point *ereg( P ), say 




! : 
i 

i : 
I ' 
i 



■I- 



^V^i i ■^■■J^ ^ 5' W. i - ' 



■. . T ■ 

' . ..[ I 



into x* - + 8{x,p)}' 



SiiTta' 4 ^ 3°1.5; F U the idenUty 

function. 



T * us study what happens in this pro- 
Let us stuu? invariant under . 
cess. Clearly, sites are in d ^ 
Farh separator sep(p, 4) 1S Qeiu . 

toaVert lS S Consequently, the 
same x 2 -coordinate^ in . 

deformed region «g bounded half 
tersection ^ upper en- 

planes of R ■ ^ re tV t * preS erves the 
velope model shows that p ^ 
topological Properties of a Voro 
gram. For instance, the interior 





AC M Costing Surveys, Vd. 23. No. 3, September 1991 



t 

e 
t 

r 
e 
f< 
n 



w 
m 
fi 
tl 

s: 
ti 
s; 
c 
ii 
a 
t; 

o 
h 
r 
c 
v 




Voronoi Diagrams 



395 




\ K 



Figure 41 . Deformmg a Voronoi diagram [Fortune 
1987]. 



Lebesgue measure zero ^ 

ad^ntage ^f™ mmost point of 
orlgl nal is that Excep . 

each re g 10 V„ n f for the boitommost sites 
n^freSart are unbounded^ 
Sire If shows the Voronoi diagram i of 
Ses (left) and its deformation (right). 



3.3.2 Plane-Sweep Technique 



The properties of V*(S) are now «- 
, -i. j L itcs construction using the 

con teucted from below by keeping ^the 
nvariant that the portion °* the <J*ect 
below L is complete at any pom i 
During the plane sweep, the cross section 

lem instead of a two-dimensional static 

^"applicability of the plane-s-eep 
tecti^crucially 

Voronoi diagram of S is to be > 
in its original: A site p e S migh , he 
above L while its region reg(p) alreaoy 
xlds to below L. And it i» «gj££ 
to know when regions begin to intersect 
7 in order to integrate them properly 
fnto th" cross section. We do -t run into 
such problems when constructing V (S) 
however, since the bottommost point of a 
deformed region is just its defining site. 



Initially when L is below S, the por- 
l r v*(^ below L contains no de- 
f n i edge Hence the cross section is 
formed edge . d the cross 

Cm Pon has to boated* either by start- 
section has to up ^ 

ing a new "^^^"a^rtei 

^fcorsTom MSecting two de- 

that comes uw nroperties of 

SB? dSSlXbovTwiKsure the 
c ^2tne e S s Cr of I^.^^SSS-?. 

Son) in a a ? d 

SSrt in a Priori queue ordered by 
? "coo d nates" There are at most n + 
/o critical points since V(fa), ana 

has at most 2 

rf V(S) twice. A. opdat^ »^ ° J th 
constant number of <V°™ K °t£ l mst 
data structures each of ""icn „ 

computed >»0(» 1 »8 t t; the origina , 
d tSStely° ffi? —ion 

translates au intersection of 

fOT wit is then a circle with 
cone(p) with « d indeed 

center p and radius _ "JP^j ag the 
this type may ^ jnterpre tea 

f rlearlv remains unatleciea. 
of cieariy i cll : tP d for a plane- 

tossing ^ 

The distance of a pom . > from » 8 
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as to have the lower endpoint of that 
segment as its bottommost pain j. Al- 
though details are more complicated a 
reasonably simple 0(nlog ») tune and 
^ n\ space algorithm is obtained. In par 
ti ular the media/ axis of a simple pdy- 
gon fs constructable via P^wee?: ^ 
Sges of the polygon are taken as line 
sefment sites, and the parts of the dia 
Sim outside of the polygon are ignored 
Ethat the separator of two segments 
consists of portions of straight lines and 
parabolas; compare Figure 25. 

3.3.3 Constrained Voronoi Diagrams 

A further type of Voronoi diagram, and 
I t which the Plane-sweep technique 
applies particularly nicely, is the c assi 
cal diagram constrained by a set ot line 
sJmenYs [Seidel 19881. : Let S denote , a 
Zt of n point sites in R and let T be a 
se of Uncrossing line segments with 
their endpoints in S. Note that m = 0(a) 
S e thfsegments form a planar ^ 
We view the segments m T as obstacles 
Id define the bounded distance between 
two points x,yeR 



as 



{ 5(x, y), 



00 



if xy n T = 0 , 
otherwise . 

The b-regions 
reg 6 (p) 

= (*£.R 2 |&(x,p)sK*.9M eS i 



\^7S^sre generally 
nSonvex (near the corresponding ob- 
Sefwhereas 6-regions of 
sites have to be convex. Both typeset 
6-regions, and necessarily the , former 
4. «<T micrht have portions ot oostacieb 
S'han o? P erpe P ndicular bisectors ^ 
edges Clearly, V,(S. T) is a polygonal 

"TheToS Voronoi diagram, is now 
SLd to the co— b Vo = 
diagram, V c (o, i ), " l ^ r the 5. 

ifying, for each obstacle teJ, the 



regions it supports Let t support , 

Z t„ \ reed pJ from the left and 
r6 i 6 S 1 ' ' ' ' Veg b S from the right. The 
forme^iregioS a?e extended to the left 
of t as if only sites Pl , . . • , P„ were pres. 
ent and the latter 6-regions are ex- 
tended to the right of Us if only sites 
a Qi ™ ere P resent - Note that Pi = q, 
has to occur twice, namely for the end- 
points r and s of t. So reg fc (r 
are extended to both sides of t. The c-re- 
Zn reg (p), of a site p e S is defined as 
fhe union of reg 6 (p) and its extensions 
riri, rwoect to all obstacles, c-regions 
d^l^mV ovilap, so they fail to define 
fcell complex and the graph , defmed by 
their edges is nonplanar See If igure -.U 
iflStrating the bounded Voronoi dia- 
Sam Solid) and the constrained Vorono, 
Sam (solid or dashed or dotted) for 11 
sites and 2 obstacles (bold) 

The significance of V e {b>, 1 ) is aue 
to 2 dual structure, the constrained 
Delaunay triangulatwn ot b ana 1. 
fntuTtive'ly, this triangulation is as sim 

lar as possible to the ^ica 1 .^ 1 *^ 
triangulation but integrates the obsta. 

sTr, a structure with -everalp^ 
cal applications. Formally, it consists 01 
he obstacle segments in T and in add, 
tion, of all edges between sites p.eS 

that have ^ ^at 

circle enclosing only sites re ^ 
least one of b(r, p), o(r, q> 

Vo^ol dSr^Note f that ^tabjdjr 
of the endpoint sites of an obstacle may 
contribute two edges (on e ther side 01 
he obstacle) rather than on lyj 
S contains n sites, T contains m obsta 
cles, and no triangulation of n points can 
have more than 3n - 6 edges, 
conclude that V e (S, T) realizes at most 
1 r, + m - 6 edges and thus has a size ui 
O(n) Figure 27 shows the constrained 
Dilaunay^riangulation that corresponds 
£ ^ ^ constrained Voronoi diagram m 

Fi Crus 2 now come to the plane-sweep 
conduction of VJLS. T). To aid the 




ACM Computing 



Surveys, Vol. 23, No. 3, September 1991 




Voronoi Diagrams 



397 




Flaure 42. Bounded and constrained Voronoi dia- 
gram [Seidel 1988]. 



points arising during the construction of 
V?{S T) and all V* is proportional to 
the size of V C (S, T) and thus is O(n). 

Using appropriate data structures sup- 
porting the m + 1 plane-sweeps allows 
us to implement the operations at each 
critical point in O(log n) time so an 
OinXogn) time and O(n) space algo- 
rithm for computing V C (S, T) « ob tamed. 
In particular, the constrained Delaunay 
triangulation of S and T is computable 
within these bounds since it can be ob- 
tained in O(n) time from V e (S, 1 ) by 
exploiting duality. 
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intuition, we think of V C (S T) as being 
embedded into m + 1 parallel p lanes in 
the following way. One plane, call it « 0 , 
contains the bounded Voronoi diagram 

V (S T) Each of the m planes, H t , tor 
the obstacle teT contains the diagram, 

V induced by the extensions with re- 
spect to t for all o-regions supported by 
t Clearly, this embedding comprises the 
same information as is inherent in 
V(S T) The idea now is to sweep « 0 
and 'all H t simultaneously by (essen- 
tially one and the same) horizontal line 
L in order to construct V 6 (S,T) and all 

V for teT, separately. Again these m 
+ 1 diagrams are not c 

rectly but via their deformations 
Y*(S T) and V* that come from map- 
ping upward, within the respective plane 
each point x inside an (extension of a) 
6-region reg b (p) by 5(x, p). Since an ex- 
tended 6-region is just a classical Voronoi 
region, its deformation has the defining 
site as the bottommost point (unless the 
site itself is bottommost for the respec- 
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Voronoi diagram, regions are created as 
sites are hit by L, and edges are created 
as L hits intersections of deformed bisec- 
tors. One must, however, take care that 
edges are assigned correctly to their em- 
bedding planes. When a deformed edge 
in H 0 meets a deformed obstacle t , then 
its construction is continued within H, 
and vice versa. The number of critical 
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me classification and prediction of the ro^s for ^"^^ 
resentation of vanous reagents independent of the ^^^^^^It a'self-organizing mode! 
the FRAU (Field-characterization for Reaction ^^^XS^U^ counter-propagation (CP) type of 
considering both negative and absent data was ■^^^ZXeno*^ CP Lhonen neural network 
=^e:« = — model *at gl ves good answers predicting 

roles of the reagents. 



Reagents play important roles in chemical reactions, and 
very various reagents have been developed and used in syn- 
thetic studies. The roles of reagents may vary according to 
the substrate, solvent, and so on. In other words, a reagent 
has various potential roles in reactions. If the roles of the 
reagents are numerically predicted on a computer before an 
experiment, it might be very useful in broad fields of syn- 
thetic study. For example, the most preferable reagents lor 
the desired reaction can be chosen before synthesis . 

In general, reagents in similar structures and electronic 
features have similar roles, because the potential roles ot 
a reagent are largely related to the features. If sufficient 
methods to numerate the structural and electronic features 
and to explain the roles based on the numeric features are 
possible, numerical prediction of the roles can be achieved. 

Recently, we developed the FRAU (Field-characteriza- 
tion for Reaction Analysis and Understanding) system and 
demonstrated its usefulness to discriminate similarities and 
differences in structures as well as the roles of metallic 
reagents based on FRAU's features for metallic atoms and 
atoms connecting to the metallic atoms. 1 ' 2 

FRAU estimates the possibilities of occurring reactions 
for a molecule (e.g., a reagent) based on electrostatic and 
steric interactions with a pseudoreactant (e.g., a substrate). 

# PRESTO, JST(Japan Science and Technology Corporation). 



The good points of FRAU are to take account of three-dimen- 
sional field around a molecule, and to measure the molecular 
features, called FRAU features (FFs), based on interactions 
with a pseudoreactant to analyze chemical reactions. The 
input data of FRAU are a molecular structure with atomic 
charges. An arbitrary computational level to obtain the in- 
put data is adaptable according to the purposes; namely, the 
degrees of accuracy and the calculation time are adjustable 
to the purposes. Thus, FRAU makes it possible to rapidly 
obtain more accurate properties for molecules than the prop- 
erties based on topological relationships between atoms. 

The number of the FFs, which are estimated for each 
atom depends on the number of atoms in a molecule. The 
atomic features are necessary to know important sites and 
directions for occurring reactions. However, in order to 
compare between reagent molecules, representation by the 
same number of features for various reagents independent 
of the number of atoms are needed. The same dimensional 
representation of molecules is one of the important problems 
in reaction classification studies, where comparing points are 
needed to discuss the similarities and differences. 

In the meanwhile, a neural network is one of appropriate 
methods for treating multi-dimensional nonlinear data. We 
used neural networks for reaction classification studies • be- 
cause the factors controlling reactions can be considered as 
data in multi-dimensional space and a reaction resulting from 
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complicated interactions among these factors can be consid- 
ered s nonlinear relationships of these mulu-dunensional 
Sa Pioneer studies that used neural networks for reasons 
w re P formed by Gasteiger et al.- The first one was an 
^plication of a multilayer back-propagation neura«k 
and an associative memory system to learning and predicting 

neural network introduced by 
Teu^ Kohonen-'^has been applied » 
tion studies, and its usefulness has been m . 
The Kohonen neural network makes .t possible to class ty in 
nut data based on the nature being inherent in the input data 
Th s i er is suitable for classifying reactions without 
^ pr conception of categories or types of — s as *e 
templates This is one of the reasons why we have been using 
rSnen neural network for our reaction classification 

^e'lssification can be developed for predictions by 
a counter-propagation type of 

fhereafter called CP Kohonen neural network), that is a 
Kohonen Neural network to which known signals and an out- 
putTyS are added. The processes of 
sification and prediction are formally similar to the thinking 
ty oi "a chemist's brain, where, following the recognition 
Zl classifications of reactions, predictions are performed 
The CP Kohonen neural network gives an organizing map 
There input data are classified in the basis of the s.milantie 
I Z nput data with a map, where the known signals are 
trl d according to the similarities in the input data. Good 
points of the CP Kohonen neural network comparing a back- 
pr P agat,on type is to give a stable interpolated answer and 
Tshow us maj>r factors of the input data that concern with 
he tested similarities because the known signals do not 
nnuence the classification. The CPKo = -^ 
work has an advantage to treat ^ul^LZ 
ourposes although there are many statistical methods used to 
Si major'factors, such as principal «**c^"£ 
vses The CP Kohonen neural network trams values for 
p diction according to similarities in the 
f the similarities in the input data have Bood^^ 
with the similarities in the known signals, the CP Kohonen 
Network produces a good model ^ ^ZZ 
nrooriate answer. However, in the onginal CP Kohonen 
network, which treats one type of knowr .signal, U 
is not distinguishable between negative and absent data in 
ord to treat reaction data, it is needed to treat absent dam 
Gently from negative ones, e. g., reaction, ^reag 
that did not occur and work, respectively. The negative 
are important to consider chemical reacuons. 

The purpose of this article is to develop and to show a 
processfoTthe characterization, classification, * md pre .tac- 
tion of the potential roles of reagents usmg comical in 
Ration w'h more accuracy than the properties ; ased on 
the topological relationships between atoms. For this pur 
nose we first solved problems about the same dimensional 
;.sIntation of reagents based on FFs 
number of atoms. Second, we improved the CP Kohonen 



neural network to treat absent data differently from negative 
ones We men have used these results to construct a reagent- 
"les c^el2n model that predicts the potential roles of the 

reagents. 



Method 

(Vi„ \) using the Kohonen neural network. Then, an origin 

mles hv the model were examined for the reagents aiiu 
of reagents as test data (Hg. 2). Detailed procedures are described 

te 2 ^odel Structures. For the first purpose that is represen- 
t L ^molecules in various structures by the same number of 

SSW, Grignard reagent, 
L 1 were used as the models. The second purpose is to predict 
S potS -its of the reagents from .he reagents' = by a 

( , g ., solvents, *^* K >™™^Z 7i 1 rive been 
"ZZ J^^^K--^ Boranes are 
i:„tKrn Senc structures, and it is coasted thatthe 

SJS-cfc, in hydrations are —^^013^) 
i, aaak such as diisobutylalununium hydnde (DIBAL) 
STe £c iSSd" bi monomenc active species in hydro- 
(15) are also consio equilibrium mixture 

aluminauons. Gngnard reagems e c i ear solvated 

which may .^-^^^• - t^e fl -^ «• 

C ^Z^ZtZ:^ because they are vaned by 

s r^s« »» ~ it,; •srr 

Kp rnnsidered to be essential to determine the roles It me sun 

next stage. The structure s w e v geometries and 

r*\ mnlf-rnlar-orbital calculations to ODiain uic B 
( ) molecular oro . ido ca i clulati0 ns were used 

charges for FRAU calculations consisted met allic atoms, 

here to obtain reliable properties for toe consist ea ™ 

3. Classification of Reagents. Calculate of FF s FRAU 
, w„t FFc for the thirty kinds of reagents. FRAU estimates 
calculated Fhs tor me inmy electrostatic 
three kinds of FFs: extent of reacuon field (r-rMcO, eicx 
* ? ) and steric feature (FF steri c). FF eleOT > and FF ;K hc 

feature ™ d S interactions with unit charge (+1) 

were calculated in the basis oi u«b« ™ n ectivelv. 

were described in our previous paper cjd _ 
Selection of Representative FRAU Futures ^ rf 

culated to each atom in a molecule. Thus, the nuro 
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Fig 1 Model reagents (planar descriptions). Three dimensional structures were optimized by ab initio RHF/3-21G( ) calculations. 
Optimized geometries for 12, 13, and 14 were octahedral type of structures. 
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Fig. 2. Reagents for test. They were not used for construction 
of the model. 

the three kinds of FFs is the same as the number of atoms in the 
molecule. For example, for dicyclohexylborane, which contains 36 
atoms (C12H23B), 36 sets of the three kinds of FFs are obtained. 

In order to compare between reagents in various structures, it 
is needed to represent the reagents by the same dimensional fac- 
tors. Thus, the same dimensional representative FFs that have good 
correspondence with similarities in the structures and roles were 
selected after investigating various combinations of FFs. 

Classification. The selected FFs were input to a Kohonen neu- 



ral network as discriminators to classify thirty kinds of the reagents 

(Fig- D- . „ , 

The Kohonen neural network outputs a planar map that is called 

as Kohonen map, where similar input data are set into the same or 
close neuron, and does not perceive boundaries between clusters of 
the input data on the Kohonen map. The boundaries were deter- 
mined by degree of gap between the weight vectors of neighboring 
neurons by a U-matrix method. 22 

4. Construction of a Reagent- Roles Correlation Model. 
Counter-Propagation (CP) Kohonen Neural Network. A 
CP Kohonen neural network is constructed from a Kohonen neu- 
ral network by the composing of output layers and known answer 
signals (Fig. 3-(l)). During the training of a CP Kohonen neural 
network, a modification of the weight vectors between the active 
and output layers is also performed. In the training, the position 
of the connection between neurons of active and output layers is 
determined according to similarities in the input signals, and the di- 
mension and value of the weight vectors between output and active 
layers are determined according to the signals of a known answer. 
Weight vectors between output and active layers corresponding to a 
known answer are trained according to the similarities in the input 
signals. 

Modification of the Original CP Kohonen Neural Network. 

In order to treat reaction data, negative and absent data should 
be distinguished; namely, it should be distinguished whether the 
reaction (reagent) did not occur (work) or the data is just absent 
from the data set. The distinction is hardly adapted to the original 
CP Kohonen neural network that has only one type of known an- 
swer signals. Thus, the original CP Kohonen neural network was 
modified according to what follows. 14 

The architecture of the modified CP Kohonen neural network is 
shown in Fig. 3-(2), where known answer signals for the presence 
of data are added. The additional signals with the other known 
answer signals control the training rate of weight vectors between 
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twork. (1) An original ^P^^^^^'rf weight vectors 
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Construction of the CF n ™ , d t0 construct a 

The modified CP Kohonen neura ""7* "reactions. In the 
correlation modelbetween reagents ^*^ m ■ t signa l s , 
construction, the selected representative FFs as the p 



first known answer signals, and 
tbe roles of the o^TL as the second known 

the information on the presence o mne 
answer signals were input As the inforrnau on on the 

types of roles of the reagents were set wim 
plence of the ^^^XXnoZo^ for each 
or 0), as shown in Table 1 AH <* ^ roles 

reagent are listed in Table 2. The input u 23.24 ^ we ii as 

determined front some « ^ g ^ 
SYNLIB* and ISIS da ^ aseS , \ binary notation was 

again st the «-^^ m nS^Tni^ necessa^ for 

determined as follows, u <m * 
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Table 1 . Binary Notations Used as Input and Known Signals in the Training of CP Kohonen Neural Network 

There is no data. 



Reduction of ketone to alcohol 
Reduction of aldehyde to alcohol 
Reduction of carboxylic acid to alcohol 
Reduction of ester to aldehyde 
Reduction of ester to alcohol 
Reduction of epoxide 
Hydroboration or hydroalumination 
Base 

Alkylation 



There is data where 
the reagent worked. 



There is data where 
the reagent did not work. 



Known 
signal 



Presence 
of data 



Known 
signal 



Presence 
of data 




Known 
signal 



Presence 
of data 



0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 



Table 2. Input Signals for Thirty Reagents in the Training of Kohonen Neural Network a ' b) 



Reagent Reduction of Reduction of 
No. ketone aldehyde 



to alcoho 



1 

2 
3 
4 
5 
6 
7 
8 
9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 



(1) (2 



0 
0 
0 
0 
0 
0 
0 
0 
0 



to alcohol 



Reduction of 
carboxylic acid 
to alcohol 



Reduction of 
ester 

to aldehyde 



Reduction of 
ester 

to alcohol 



Reduction of Hydroboration or Base Alkylation 
epoxide hydroalumination 



0 
0 
0 
1 
1 
0 
1 
1 
1 



0 

1 
1 
1 

1 

0 
0 

1 

0 

1 
1 
1 
1 
1 
1 

0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 



0 

1 
1 
1 
1 

0 
0 

1 

0 

1 
1 
1 
1 
1 
1 

0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 



0 
0 

1 

1 

0 
0 
0 
0 
0 

1 

0 
0 

1 

0 

1 

0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 



0 

1 
1 
1 

0 
0 
0 
0 
0 

1 

0 
0 

1 
1 
1 

0 
0 
0 
0 
0 
0 
0 
0 
0' 
0 
0 
0 
0 
0 
0 



0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 

1 

0 

1 

0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 



0 

1 

0 
0 
0 
0 
0 

1 

0 

1 

1 

0 
0 
0 

1 

0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
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0 
0 
0 
0 
0 
0 
0 
0 

1 

0 
0 

1 

1 
1 

0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 



0 

1 

0 
0 
0 
0 
0 

1 

0 



0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 



0 
0 
0 
0 
0 
0 
0 
0 
0 

1 

0 
0 

1 

0 

1 

0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 



0 
0 
0 
0 
0 
0 
0 
0 
0 

1 

0 
0 

1 

1 
1 

0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 



(1) (2) (1) (2) (1) (2) (1) (2) (1) (2) (1) (2) 



1 
1 
1 
1 

0 

1 

0 

1 

0 

1 

0 
0 

1 

0 

1 

0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
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0 
0 
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1 
1 
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0 
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(1) (2) (1) (2) 



0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 1 

0 0 0 

0 0 0 

0 0 0 

1 1 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 

0 0 0 

1 1 



0 
1 0 
1 0 



0 0 1 

0 0 1 

0 0 1 

0 0 1 

0 0 1 

0 0 1 



0 
0 
0 
0 
0 
0 
0 
0 

1 

0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 



a) Reagents' numbers that correspond to those in Fig. 1. b) Roles in reactions. (1) Known answer 
answer signals for presence of data. 



signals for the roles of reagents. (2) Known 



an interested reagent to work, such as the reduction of epoxides 
under NaBFLi with BF 3 *Et 2 0, it was not determined whether the 
reagent works, because in such a case the interested reagents are 
changed to the actual active species. If an interested reagent works, 
even against a specific substrate, it was determined that the reagent 
works. The volume of the voice for absent data was set at 30% of 



that for the existing data. 

5. Prediction of Roles of Reagents. Reproduction and pre- 
diction by the constructed model were investigated for ten reagents 
used in the modeling (structures of 4, 6, 8, 10, 13, 15, 16, 21, 
23, and 27), and for ten reagents not used in the modeling as test 
data (Fig. 2), respectively. Reagents 38, 39, and 40 have not been 
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reported, and for the others, experiments have been reported. 



Results and Discussion 
1 Results of the Classification. A set of six parame- 
ters consisting of maximum and minimum FF e .ectro values in 
a molecule with values of FF field and FF steric on the sl tes hav- 
ing the maximum and minimum FF elect ro successfully distin- 
guished similarities and differences in the reagents and was 
SLed as a representative set of FFs. The selected repre- 
sentative sites and the FFs for each molecule are shown » 
Fie 5 and Table 3, respectively. 

The selected set of FFs was used as discriminators to clas- 
sify the reagents by the Kohonen neural network. The ; rest- 
ing Kohonen map is shown in Fig. 6-0). In F>E- 6 ;®' ^ cn 
quale is a neuron. Reagents were set on neurons where they 
were drawn. The black solid lines are boundaries recognized 
by he U-atnx method. The bold solid lines mean a tager 
difference between neurons than the fine sohd lines, and the 
fine solid lines mean a larger difference between the neurons 
than the dotted lines. The actual shape of this map is a toms. 
For visualization, the torus is cut along two perpendicular 
Unes and the surface is spread into a plane. Thus, the top of 
Ais map connects with the bottom, and the left edge connects 

with the right one. • m ;wFF<5 are set 

On the Kohonen map, input data having similar FFs are set 
into the same or close neuron. Thus, the correlation between 
the FFs and types of the reagents can be found by labeling 
the map according to the similarities in the types Figure 6 
is colored according to the similarities in the roles of the 
reagents It shows that reagents which play similar roles 



were successfully set on the same or close neurons forming 
clusters. Borane-dimethyl sulfide (BMS) (10) was isolated 
from cluster b, colored grayish blue. The location of BMS 
is suitable because the roles of BMS are rather similar to 
those of boranes (cluster a colored light blue), even though 
BMS belongs to borane complexes. The results show tha 
the selected set of FFs has a good correspondence with the 
similarities of the roles of the reagents. 
. The reason why several combinations of maximum and 
minimum values of FFs were mainly investiga ^ .eta* 
ing representative FFs is because it ,s considered mat to.se 
xLme values might have a hxgher influence to detennine 
the characters of the reagents. A good correspondence of ±e 
six parameters (the maximum and rmnimum FF dectro value 
in a molecule and values of FF fieM and FF steric on the si e 
having the maximum and minimum FF electr o) with the roles 
shows that those sites having the maximum and minimum 
electrostatic interaction energy are major for determining he 
reagents' characters. Then, in order to know which .are -.toe 
Jor parameters to discriminate the simtlanUes and differ 
ences in the reagents, neurons are colored according to the 
size of weight vectors corresponding to the six Parameter^ 
as shown in Fig. 7. Referring to the maps .n Rg. 7 w.ft 
the classification results in Fig. 6 let us know that the ma- 
or factors to detennine the boundaries are the maximum 
ind minimum of FF^o values. Various combination .of 
FFs were investigated to select the representative FFs, and 
a set of the maximum and minimum FF electro values was no 
sufficient to clearly discriminate the similarities, although it 
wt • To investigated. This means that the steric interaction 
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5 10 15 , • • pf values were important sites to determine the 

Flg 5. Selected representative sUes. The sites having maximum and ~m ^ ^ rf ^ 

8 re gents' characters. A set of six parameter, of maximum and that distinguishes similanttes and 

FfL on the sites having the maximum and minimum FF Btaaro were selected 
differences in the reagents' roles. 
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Table 3. 



Selected Representative FRAU's Features (FFs) for the Reagents 



Reagent 
No. 



Sites having minimum FFeiectro 



Sites having maximum FFeiectro 



FF 



electro 



FF ; 



steric 



1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 



-13.454 
-9.1724 
-5.2507 
-11.480 
-8.1680 
-11.526 
-33.422 
-26.279 
-35.407 
-19.256 
-39.802 
-67.907 
-29.474 
-55.343 
-20.992 
-40.573 
-40.753 
-36.978 
-29.897 
-35.948 
-36.011 
-25.630 
-25.464 
-22.381 
-23.125 
-24.032 
-24.554 
-25.820 
-23.023 
-20.703 



256.961 

271.503 

267.347 

238.851 

190.646 

235.594 

222.098 

215.623 

430.112 

260.357 

269.070 

275.703 

239.217 

261.024 

241.131 

231.209 

230.966 

230.305 

230.305 

230.312 

230.151 

209.678 

209.524 

208.736 

208.351 

208.623 

208.811 

209.217 

208.623 

208.659 



FFfield 

253068 
252849 
248989 
251560 
524012 
251559 
284007 
326616 
68969 
250266 
256267 
370214 
257019 
253623 
328616 
362756 
362876 
362680 
360838 
362154 
361815 
946002 
945928 
945287 
946616 
946875 
946937 
946016 
946772 
945317 



a) Reagents' numbers that correspond to those in Fig. 1. 

energy and field extension around the sites are also impor- 
tant to determine characters of the reagents as well as the 
electrostatic interaction energies on the sites. 

2. Results from the Reproduction and Prediction. 
Reagents Used in the Modeling. The predicted degrees 
of the roles for ten reagents used in the training of the CP 
Kohonen neural network model are given in Table 4 along 
with the input binary notations that stand for the roles and 
presence of the data. In Table 4, a value that is described 
on values in parenthesis is the predicted one; the first and 
second values in parenthesis are the input signal for the roles 
and for presence of the data, respectively. The range of pre- 
dicted values was from 0.0 to 1.0. The values of 0.0 and 1.0 
correspond to the lowest and highest possibilities to react, 
respectively. 

Almost all of the input values were successfully repro- 
duced, as shown in Table 4, where the predicted values that 
are largely different from the input ones are underlined. This 
shows the good capability of the model to reproduce the input 
values. The underlined cases are discussed below. 

The predicted values for roles of isopinocamphan-J-y - 
borane (6) as a reducing agent from aldehydes to alcohols 
and as that from carboxylic acids to alcohols were extremely 
different from the input signal. The input signal for the role 



FFeiectro 

10.919 
12.776 
13.274 
14.973 
44.590 
7.5380 
15.621 
16.355 
24.953 
14.445 
35.413 
26.285 
41.683 
35.444 
29.007 
26.975 
25.117 
33.482 
42.365 
32.833 
32.636 
38.978 
38.632 
46.611 
50.338 
49.068 
48.018 
42.793 
53.404 
51.967 



FF 



stenc 



FFfield 



281.860 
292,972 
313.062 
477.965 
203.644 
424.580 
285.051 
261.958 
304.273 
262.250 
258.575 
19.6603 
99.0511 
39.2887 
265.351 
250.520 
250.187 
228.767 
211.166 
211.975 
211.224 
271.667 
271.077 
244.503 
224.836 
236.137 
246.892 
241.687 
237.187 
231.063 



309377 
292352 
151018 
18533 
438130 
48907 
216902 
238697 
206290 
233260 
258254 
2422772 
847636 
1526515 
552630 
678728 
704083 
779961 
773137 
772059 
771967 
594481 - 
615552 
704713 
703261 
680578 
670890 
642569 
692784 
698610 



of 6 was 0 and the predicted value was 0.8. Namely, the 
organizing map predicted 80% of possibility for 6 to work 
as reducing agents from aldehydes to alcohols and from car- 
boxylic acids to alcohols, although the input signal denoted 
that there were no data. The predicted degrees of the roles 
were according to the similarities in the input FFs between 6 
and the other reagents, such as 1,1,2-tnmethylpropylborane 
(thexylborane) (4), gathered closely in the Kohonen map 
The predicted values for 6- are suitable to knowledge that 

chemists possess. . 

For LiAlH 4 (13), the predicted possibility as a reducing 
agent from esters to aldehydes was 0.4, although the input 
signal denoted that it does not work. This was caused by 
nearby position of 13 to that of 15. Incidentally, it is known 
that 13 with Et 2 NH reduces esters to aldehydes. 

For DIBAL (15), 1.0 was predicted to reduce carboxylic 
acid as the possibility. It is known that 15 does not reduce 
carboxylic acids, except pyrazolecarboxyhc acids. In this 
execution, for a reagent that works even against » a specific 
substrate, it was determined the reagent works. Thus, this 
result is an instance of an overestimate. Some extension is 
needed in the future, where information on the majorities ot 
the negative data, namely, reaction data on 15 not reactive 
to many carboxylic acids, are also considered in training the 
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511 £a| boranes 

[53 borane complexes 
hydrides 

^ reducing agents containing N atom 
a reducing agents containing O atom 



bases 

Grignard reagents 



______§ _ J Each small square in the map 

R^^^^ SS'SS^^ — ■ B,ack solid lines r.T^Z 

is a neuron. Neurons obtaining the input data are ^ to f^^ e neurons than the fine solid lines, and the fine sol d 
of 31 _40 were mapped on the neurons pointed by arrows. 



'electro 



For sites where 
minimum FF^ectro 
values were obtained 




FFsteric 





For sites where 
maximum FF e iectro 
values were obtained 



BBHWWWwKSig^ l parameters of FFs. Correspondence 

Fig. 7. ^—-^X^ 



values. 

model A predicted possibility as a base for 15 was 0.5 but 
r h einput signal denoted there was no data. Thts ts suttable 
to the nature of 15, which has the basisity. 

"edicted value for the roles of neopentylmagnestum 



chloride (27) as base was rajer different _from * eu, u 
signal. The input signal for the role for 27 was 0 _an I m 
predicted value was 0.6. The predicted degree of . the ^ roles 
was according to the similarities in the mput FFs between 27 
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Table 4. Results from Prediction of Roles of Reagents Used in Construction of the Model 

Reduction of Reduction of Reduction of Reduction of Reduction of Reduction of Hydroboration 
ketone aldehyde carboxylic acid ester ester epoxide ornydro- 
to alcohol to alcohol to alcohol to aldehyde to alcohol alurrunation 



Base Alkylation 




8 



H / 



10 



LiAIH4 



13 



15 

Mg 

* 16 

21 

Cl 23 



27 



1.0 

(1,1) 


0.8 

(1,1) 


0.8 
(1,1) 


0.1 
■ (0,0) 


0.1 
(0,0) 


0.1 

(0,0) 


1.0 

(1,1) 


0.0 

(0,0) 


0.0 
(0,0) 


1.0 

(1,1) 


0.8 
(0,0) 


0.8 
(0,0) 


0.1 
(0,0) 


0.1 
(0,0) 


0.1 
(0,0) 


1.0 

(1,1) 


0.0 
(0,0) 


0.0 
(0,0) 


1.0 

(1,1) 


0.9 
(1,1) 


0.1 

(0,0) 


0.0 
(0,1) 


0.0 
(0,1) 


0.1 
(0,0) 


0.9 

(1,1) 


0.1 

(0,0) 


0.0 

(0,0) 


1.0 

(1,1) 


0.9 
(1,1) 


0.7 
(1,1) 


0.0 
(0,1) 


0.7 

(1,1) 


0.7 

(1,1) 


1.0 

(1,1) 


0.0 

(0,0) 


0.0 
(0,0) 


1.0 

(1,1) 


1.0 
(1,1) 


0.8 

0,1) 


0.4 
(0,1) 


1.0 

(1,1) 


0.8 

(1,1) 


0.9 

(1,1) 


0.8 

(1,1) 


0.0 
(0,0) 


1.0 

(1,1) 


1.0 

(1,1) 


1.0 

(1,1) 


0.8 

(1,1) 


0.9 

(1,1) 


0.9 

(1,1) 


0.9 

(1,1) 


0.5 
(0,0) 


0.0 
(0,0) 


1.0 

(1,1) 


0.0 

(0,0) 


0.0 
(0,0) 


0.0 
(0,0) 


0.0 
(0,0) 


0.0 
(0,0) 


0.0 
(0,1) 


0.0 
(0,0) 


0.0 
(0,0) 


1.0 

(1,1) 


0.1 

(0,0) 


0.1 
(0,0) 


0.1 

(0,0) 


0.1 
(0,0) 


0.1 

(0,0) 


0.0 

(0,0) 


0.0 
(0,0) 


0.0 
(0,0) 


0.1 

(0,0) 


0.2 
(0,0) 


0.0 

(0,0) 


0.0 
(0,0) 


0.0 
(0,0) 


0.0 

(0,0) 


0.0 
(0,1) 


0.9 

(1,1) 


0.4 

(0,0) 


0.0 
(0,0) 


0.0 
(0,0) 


0.0 
(0,0) 


0.0 
(0,0) 


0.0 
(0,0) 


0.0 
(0,0) 


0.0 
(0,0) 


0.6 
(0,0) 


0.9 

(1,1) 



and the other base reagents gathered closely in the Kohonen 
map. This result is also suitable to the nature of 27, which 
has the strong basisity, although it is not actually usually used 
as base. 

These results demonstrate that the degrees of similarities 
in the input FFs have a good correspondence with the degree 
of similarities in the roles of the reagents, and that the CP 
Kohonen neural network model was successfully trained. 

Reagents Not Used in the Modeling. In the CP Kohonen 
neural network model, ten reagents not used in the modeling 
(Fig. 2) were allocated as shown in Fig. 6-(II), and the pre- 
dicted degrees of the roles are listed in Table 5. In the first 
column of Table 5, the reagent structures and the structural 
nos. with q-error values in parenthesis are shown. The q- 
error values indicate a gap between a reagent's FFs and the 
weight vectors of a neuron where the reagent was mapped. 
Thus, lower values of the q-error denote a higher reliability 
of the predicted values. 

In Fig. 6, reagents were allocated to the neurons indi- 
cated by arrows. The reported borane of 9-borabicyclo- 
[3.3.1]nonane (9-BBN) (37) and ammonia- borane (36) are 
allocated into groups of boranes and borane complexes, re- 



spectively. Unreported boranes 38, 39, and 40 are allo- 
cated into groups of boranes. LiBH 4 (35) was allocated to 
a neuron neighbor to L1AIH4. MeOMgH (32), 2,6-Me 2 -c- 
C 5 H 9 NMgH (31), «-Pr 2 NMgCl (33), and MeMgCl (34) were 
allocated to groups of alkoxymagnesium hydrides, dialkyl- 
aminomagnesium hydrides, dialkylaminomagnesium chlo- 
rides, and Grignard reagents, respectively. The results show 
that each of ten reagents was successfully mapped into a 

suitable neuron. 

As shown in Table 5, for 31 and 32, the possibilities of 
reducing agents from ketone to alcohol were predicted as 
1.0. They were reported as reducing agents from ketone to 
alcohol. 27 34, which is well known as a Grignard reagent, the 
possibilities of a reagent for alkylation was predicted to be 
1.0. Low values of the q-error for them are consistent with 
these correct results. 

For L1BH4 (35), all of the predicted values, except that 
for reducing carboxylic acids, are suitable. The predicted 
values for reducing carboxylic acids as 0.8, underlined in 
Table 5, is not good, because it is known that 35 does not 
reduce carboxylic acids. For 9-BBN (37), all of the predicted 
values, except that for reducing carboxylic acids and esters, 
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Table 5. Results from Prediction of Roles of Reagents Not Used in Construction of the Model 



Reduction of Reduction of Reduction of Reduction of Reduction of of *™ ^ 

ketone aldehyde carboxylic acid ester ester epoxide or ny 

to alcohol to alcohol to alcohol to aldehyde to alcohol _ 



r 

31 (0. H 1369) 

32 (0.1242) 



1.0 



0.0 



0.0 



0.0 



1.0 



0.1 



0.1 



0.1 




0.1 



0.2 



0.0 



0.0 



34(0.1748) 

L1BH4 
35(1.0442) 

H H 

36 (0.7299) 



0.0 



1.0 



1.0 



0.0 



1.0 



1.0 



0.0 



0.8 



0.0 



0.0 



0.0 



0.0 




1.0 



0.9 



0.7. 



0.0 



1.0 



0.8 



0.5 



0.0 



1.0 



0.8 



0.5 



0.0 



39 (0.3497) 

H 




1.0 



0.8 



0.8 



0.1 



40(1.0548) 



are suitable. The predicted values for reducing carboxylic 
acids and esters as 0.7, underlined in Table 5, are not good, 
because it is known that 37 does not reduce them. The q- 
error values for 35 and 37 were higher. The q-error values 
denote a lower reliability of the predicted values. 

These results show that the predicted values are almost 
suitable to the potential of the roles, and that the reliability of 
the predicted values is given as q-error values. Thus, the CP 
Kohonen neural network model constructed here can predict 
the roles of those reagents not used in the modeling. 

Stability of the CP Kohonen Neural Network Model. 
Training of the CP Kohonen network model was repeated 
to the know stability of the model by an examination of the 
weight vectors between active and output layers. During 
the training, the weight vectors from about 0.3 to 0.7 were 
rather unstable, while those from about 0.0 to 0.3 and from 
about 0.7 to 1.0 were rather stable. The results mean that the 
predicted values by the model can be used as three or four 
classes of possibility measures. 

A good point of the CP Kohonen neural network is to know 



0.0 



0.0 



0.0 



0.0 . 0.0 



0.1 



0.1 



0.0 



0.0 0.0 



0.0 



0.0 



0.0 



0.9 0.4 



0.0 



1.0 



0.0 



0.0 



0.8 



0.0 



0.0 



0.9 



0.0 



0.3 1.0 



0.9 0.0 



0.0 0.1 



0.7 



0.7 



1.0 



0.0 0.0 



0.4 



0.6 



1.0 



0.0 0.0 



0.4 



0.6 



1.0 



0.0 0.0 



0.1 



0.1 



1.0 



0.0 0.0 



the major factors related to the roles of reagents. This is a 
different good point from a back- propagation (BP) neural 
network, and the point is the reason to use the CP neural 
network, here. A weak point of the CP compared with BP 
neural networks is that the model can not predict the roles for 
a reagent that is not similar to reagents used in the modeling. 
Thus, a combination of the good points of CP and BP neural 
networks will produce a better model for prediction, which 
is what we are planing to do. 

Conclusion 

The same dimensional representative features that were 
calculated by a FRAU system represented similarities in the 
roles of reagents in various structures. The representative 
FRAU's features shown that sites where the maximum and 
minimum electrostatic interaction energies were calculated 
are important to determine the reagents' characters. The 
roles of the reagents were successfully classified in the basis 
of the representative FRAU's features by a Kohonen neu- 
ral network. The results were applied to a neural network 
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model construction for predicting the potential roles of the 
reagents by a CP Kohonen neural network that has been 
modified to distinguish between negative and absent data. 
Thus, a reagent-roles correlation model that can predict the 
possibilities of the reagents' roles was constructed. 

The FRAU features represent more detailed similarities 
and differences in molecules and the FRAU features as dis- 
criminators for close structures and roles are now in progress 
and will be reported in the future. Furthermore, verifications 
of the predicted values by experiments are also in progress 
and will be described elsewhere. 

The benefits of a CP Kohonen neural network is to treat 
non-linear data as well as to know factors that principally 
contribute to the classification. The CP Kohonen neural net- 
work produced a reagent-roles correlation model that gives 
good reproduced and stable interpolated answers. The pre- 
dicted values can be evaluated in three to four classes. For 
more accuracy, an improvement of neural network methods 
will be needed. The reliability of extrapolated answers and 
the construction of an extended model will be investigated 
by the CP neural network model with other methods, such as 
a back-propagation neural network. 

The procedures for constructing the reagent-roles correla- 
tion model presented here are planned to apply to our reaction 
prediction studies in the future. 

We thank Dr. Kazuo Nagasawa of Synthetic Organic 
Chemistry Laboratory in RIKEN (The Institute of Physical 
and Chemical Research) for useful discussions and advice 
on the chemical reactivity of reagents. 
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